Methods for multiplex pcr amplification of target loci in a nucleic acid sample

ABSTRACT

Methods for non-invasive prenatal paternity testing are disclosed herein. The method uses genetic measurements made on plasma taken from a pregnant mother, along with genetic measurements of the alleged father, and genetic measurements of the mother, to determine whether or not the alleged father is the biological father of the fetus. This is accomplished by way of an informatics based method that can compare the genetic fingerprint of the fetal DNA found in maternal plasma to the genetic fingerprint of the alleged father.

RELATED APPLICATIONS

This application is a continuation of U.S. Utility application Ser. No. 17/018,831, filed Sep. 11, 2020. U.S. Utility application Ser. No. 17/018,831 is a continuation of U.S. Utility application Ser. No. 16/287,662, filed Feb. 27, 2019, now U.S. Pat. No. 10,774,380. U.S. Utility application Ser. No. 16/287,662, now U.S. Pat. No. 10,774,380, is a continuation of U.S. Utility application Ser. No. 16/012,667 filed Jun. 19, 2018. U.S. Utility application Ser. No. 16/012,667 is a continuation of U.S. Utility application Ser. No. 13/335,043 filed Dec. 22, 2011, now U.S. Pat. No. 10,113,196. U.S. Utility application Ser. No. 13/335,043, now U.S. Pat. No. 10,113,196, claims the benefit of U.S. Provisional Application Ser. No. 61/426,208, filed Dec. 22, 2010, and is a continuation-in-part of U.S. Utility application Ser. No. 13/300,235, filed Nov. 18, 2011, now U.S. Pat. No. 10,017,812. U.S. Utility application Ser. No. 13/300,235, now U.S. Pat. No. 10,017,812, claims the benefit of U.S. Provisional Application Ser. No. 61/571,248, filed Jun. 23, 2011; U.S. Provisional Application Ser. No. 61/542,508, filed Oct. 3, 2011; and is a continuation-in-part of U.S. Utility application Ser. No. 13/110,685, filed May 18, 2011, now U.S. Pat. No. 8,825,412. U.S. Utility application Ser. No. 13/110,685, now U.S. Pat. No. 8,825,412, claims the benefit of U.S. Provisional Application Ser. No. 61/395,850, filed May 18, 2010; U.S. Provisional Application Ser. No. 61/398,159, filed Jun. 21, 2010; U.S. Provisional Application Ser. No. 61/462,972, filed Feb. 9, 2011; U.S. Provisional Application Ser. No. 61/448,547, filed Mar. 2, 2011; and U.S. Provisional Application Ser. No. 61/516,996, filed Apr. 12, 2011, and the entirety of all these applications are hereby incorporated herein by reference for the teachings therein.

FIELD

The present disclosure relates generally to methods for non-invasive prenatal paternity testing.

BACKGROUND

Unclear parentage is a significant problem, and estimates range between 4% and 10% of children who believe their biological father to be a man who is not their actual biological father. In cases where a woman is pregnant, but relevant individuals are not sure who the biological father is, there are several options to determine the correct biological father of the fetus. One method is to wait until birth, and conduct genetic fingerprinting on the child and compare the genetic fingerprint of the child's genome with that of the suspected fathers. However, the mother often wishes to know the identity of the biological father of her fetus prenatally. Another method is to perform chorionic villus sampling in the first trimester or amniocentesis in the second trimester, and use the genetic material retrieved to conduct genetic fingerprinting prenatally. However, these methods are invasive, and carry a significant risk of miscarriage.

It has recently been discovered that fetal cell-free DNA (cfDNA) and intact fetal cells can enter maternal blood circulation. Consequently, analysis of this fetal genetic material can allow early Non-Invasive Prenatal Genetic Diagnosis (NIPGD or NPD). A key challenge in performing NIPGD on fetal cells is the task of identifying and extracting fetal cells or nucleic acids from the mother's blood. The fetal cell concentration in maternal blood depends on the stage of pregnancy and the condition of the fetus, but estimates range from one to forty fetal cells in every milliliter of maternal blood, or less than one fetal cell per 100,000 maternal nucleated cells. Current techniques are able to isolate small quantities of fetal cells from the mother's blood, although it is difficult to enrich the fetal cells to purity in any quantity. The most effective technique in this context involves the use of monoclonal antibodies, but other techniques used to isolate fetal cells include density centrifugation, selective lysis of adult erythrocytes, and FACS. A key challenge is performing NIPGD on fetal cfDNA is that it is typically mixed with maternal cfDNA, and thus the analysis of the cfDNA is hindered by the need to account for the maternal genotypic signal. Fetal DNA analysis has been demonstrated using PCR amplification using primers that are designed to hybridize to sequences that are specific to the paternally inherited genes. These sources of fetal genetic material open the door to non-invasive prenatal diagnostic techniques.

Once the fetal DNA has been isolated, either pure or in a mixture, it may be amplified. There are a number of methods available for whole genome amplification (WGA): ligation-mediated PCR (LM-PCR), degenerate oligonucleotide primer PCR (DOP-PCR), and multiple displacement amplification (MDA). There are a number of methods available for targeted amplification including PCR, and circularizing probes such as MOLECULAR INVERSION PROBES (MIPs), and PADLOCK probes. There are other methods that may be used for preferentially enrich fetal DNA such as size separation and hybrid capture probes.

There are numerous difficulties in using DNA amplification in these contexts. Amplification of single-cell DNA, DNA from a small number of cells, or from smaller amounts of DNA, by PCR can fail completely. This is often due to contamination of the DNA, the loss of the cell, its DNA, or accessibility of the DNA during the amplification reaction. Other sources of error that may arise in measuring the fetal DNA by amplification and microarray analysis include transcription errors introduced by the DNA polymerase where a particular nucleotide is incorrectly copied during PCR, and microarray reading errors due to imperfect hybridization on the array. Another problem is allele drop-out (ADO) defined as the failure to amplify one of the two alleles in a heterozygous cell.

Many techniques exist which provide genotyping data. Some examples include the following. TAQMAN is a unique genotyping technology produced and distributed by LIFE TECHNOLOGY. TAQMAN uses polymerase chain reaction (PCR) to amplify sequences of interest. AFFYMETRIX's 500K ARRAYS and ILLUMINA's INFINIUM system are genotyping arrays that detect for the presence of specific sequences of DNA at a large number of locations simultaneously. ILLUMINA's HISEQ and MISEQ, and LIFE TECHNOLOGY's ION TORRENT and SOLID platform allow the direct sequencing of a large number of individual DNA sequences.

SUMMARY

Disclosed herein are methods for determining the paternity of a gestating fetus in a non-invasive manner. According to aspects illustrated herein, in an embodiment, a method for establishing whether an alleged father is the biological father of a fetus that is gestating in a pregnant mother includes obtaining genetic material from the alleged father, obtaining a blood sample from the pregnant mother, making genotypic measurements, at a plurality of polymorphic loci, on the genetic material from the alleged father, obtaining genotypic measurements, at the plurality of polymorphic loci, from the genetic material from the pregnant mother, making genotypic measurements on a mixed sample of DNA originating from the blood sample from the pregnant mother, where the mixed sample of DNA comprises fetal DNA and maternal DNA, determining, on a computer, the probability that the alleged father is the biological father of the fetus gestating in the pregnant mother using the genotypic measurements made from the DNA from the alleged father, the genotypic measurements obtained from the pregnant mother, and the genotypic measurements made on the mixed sample of DNA, and establishing whether the alleged father is the biological father of the fetus using the determined probability that the alleged father is the biological father of the fetus.

In an embodiment, the polymorphic loci comprise single nucleotide polymorphisms. In an embodiment, the mixed sample of DNA comprises DNA that was from free floating DNA in a plasma fraction of the blood sample from the pregnant mother. In an embodiment, the mixed sample of DNA comprises maternal whole blood or a fraction of maternal blood containing nucleated cells. In an embodiment, the fraction of maternal blood containing nucleated cells has been enriched for cells of fetal origin.

In an embodiment, the determination of whether the alleged father is the biological father includes calculating a test statistic for the alleged father and the fetus, wherein the test statistic indicates a degree of genetic similarity between the alleged father and the fetus, and wherein the test statistic is based on the genotypic measurements made from DNA from the alleged father, the genotypic measurements made from the mixed sample of DNA, and the genotypic measurements obtained from DNA from the pregnant mother, calculating a distribution of a test statistic for a plurality of individuals who are genetically unrelated to the fetus, where each calculated test statistic indicates a degree of genetic similarity between an unrelated individual from the plurality of individuals who are unrelated to the fetus and the fetus, wherein the test statistic is based on genotypic measurements made from DNA from the unrelated individual, the genotypic measurements made from the mixed sample of DNA, and genotypic measurements obtained from DNA from the pregnant mother, calculating a probability that the test statistic calculated for the alleged father and the fetus is part of the distribution of the test statistic calculated for the plurality of unrelated individuals and the fetus, and determining the probability that the alleged father is the biological father of the fetus using the probability that the test statistic calculated for the alleged father is part of the distribution of the test statistic calculated for the plurality of unrelated individuals and the fetus. In an embodiment, establishing whether an alleged father is the biological father of the fetus also includes establishing that the alleged father is the biological father of the fetus by rejecting a hypothesis that the alleged father is unrelated to the fetus if the probability that the alleged father is the biological father of the fetus is above an upper threshold, or establishing that the alleged father is not the biological father of the fetus by not rejecting a hypothesis that the alleged father is unrelated to the fetus if the probability that the alleged father is the biological father of the fetus is below a lower threshold, or not establishing whether an alleged father is the biological father of the fetus if the likelihood is between the lower threshold and the upper threshold, or if the likelihood was not determined with sufficiently high confidence.

In an embodiment, determining the probability that the alleged father is the biological father of the fetus includes obtaining population frequencies of alleles for each locus in the plurality of polymorphic loci, creating a partition of possible fractions of fetal DNA in the mixed sample of DNA that range from a lower limit of fetal fraction to an upper limit of fetal fraction, calculating a probability that the alleged father is the biological father of the fetus given the genotypic measurements obtained from DNA from the mother, the genotypic measurements made from DNA from the alleged father, the genotypic measurements made from the mixed sample of DNA, for each of the possible fetal fractions in the partition, determining the probability that the alleged father is the biological father of the fetus by combining the calculated probabilities that the alleged father is the biological father of the fetus for each of the possible fetal fractions in the partition, calculating a probability that the alleged father is not the biological father of the fetus given the genotypic measurements made from DNA from the mother, the genotypic measurements made from the mixed sample of DNA, the obtained allele population frequencies; for each of the possible fetal fractions in the partition, and determining the probability that the alleged father is not the biological father of the fetus by combining the calculated probabilities that the alleged father is not the biological father of the fetus for each of the possible fetal fractions in the partition.

In an embodiment, calculating the probability that the alleged father is the biological father of the fetus and calculating the probability that the alleged father is not the biological father of the fetus may also include calculating, for each of the plurality of polymorphic loci, a likelihood of observed sequence data at a particular locus using a platform response model, one or a plurality of fractions in the possible fetal fractions partition, a plurality of allele ratios for the mother, a plurality of allele ratios for the alleged father, and a plurality of allele ratios for the fetus, calculating a likelihood that the alleged father is the biological father by combining the likelihood of the observed sequence data at each polymorphic locus over all fetal fractions in the partition, over the mother allele ratios in the set of polymorphic loci, over the alleged father allele ratios in the set of polymorphic loci, and over the fetal allele rations in the set of polymorphic loci, calculating a likelihood that the alleged father is not the biological father by combining the likelihood of the observed sequence data at each polymorphic locus over all fetal fractions in the partition, over the mother allele ratios in the set of polymorphic loci, over population frequencies for the set of polymorphic loci, and over the fetal allele ratios in the set of polymorphic loci, calculating a probability that the alleged father is the biological father based on the likelihood that the alleged father is the biological father, and calculating a probability that the alleged father is not the biological father based on the likelihood that the alleged father is not the biological father.

In an embodiment, calculating the probability that the alleged father is the biological father based on the likelihood that the alleged father is the biological father is performed using a maximum likelihood estimation, or a maximum a posteriori technique. In an embodiment, establishing whether an alleged father is the biological father of a fetus may also include establishing that the alleged father is the biological father if the calculated probability that the alleged father is the biological father of the fetus is significantly greater than the calculated probability that the alleged father is not the biological father, or establishing that the alleged father is not the biological father of the fetus if the calculated probability that the alleged father is the biological father is significantly greater than the calculated probability that the alleged father is not the biological father. In an embodiment, the polymorphic loci correspond to chromosomes that have a high likelihood of being disomic.

In an embodiment, the partition of possible fractions of fetal DNA contains only one fetal fraction, and where the fetal fraction is determined by a technique taken from the list consisting of quantitative PCR, digital PCR, targeted PCR, circularizing probes, other methods of DNA amplification, capture by hybridization probes, other methods of preferential enrichment, SNP microarrays, DNA microarrays, sequencing, other techniques for measuring polymorphic alleles, other techniques for measuring non-polymorphic alleles, measuring polymorphic alleles that are present in the genome of the father but not present in the genome of the mother, measuring non-polymorphic alleles that are present in the genome of the father but not present in the genome of the mother, measuring alleles that are specific to the Y-chromosome, comparing the measured amount of paternally inherited alleles to the measured amount of maternally inherited alleles, maximum likelihood estimates, maximum a posteriori techniques, and combinations thereof. In an embodiment, the method of claim 1, wherein the partition of possible fetal fractions contains only one fetal fraction, and where the fetal fraction is determined using the method of claim 26.

In an embodiment, the alleged father's genetic material is obtained from tissue selected from the group consisting of: blood, somatic tissue, sperm, hair, buccal sample, skin, other forensic samples, and combinations thereof. In an embodiment, a confidence is computed for the established determination of whether the alleged father is the biological father of the fetus. In an embodiment, the fraction of fetal DNA in the mixed sample of DNA has been enriched using a method selected from the group consisting of: size selection, universal ligation mediated PCR, PCR with short extension times, other methods of enrichment, and combinations thereof.

In an embodiment, obtaining genotypic measurements from the genetic material of pregnant mother may include making genotypic measurements on a sample of genetic material from the pregnant mother that consists essentially of maternal genetic material. In an embodiment, obtaining genotypic measurements from the genetic material from the pregnant mother may include inferring which genotypic measurements from the genotypic measurements made on the mixed sample of DNA are likely attributable to genetic material from the pregnant mother, and using those genotypic measurements that were inferred to be attributable to genetic material from the mother as the obtained genotypic measurements. In an embodiment, the method may also include making a clinical decision based on the established paternity determination. In an embodiment, the clinical decision is to terminate a pregnancy.

In an embodiment, making genotypic measurements may be done by measuring genetic material using a technique or technology selected from the group consisting of padlock probes, molecular inversion probes, other circularizing probes, genotyping microarrays, SNP genotyping assays, chip based microarrays, bead based microarrays, other SNP microarrays, other genotyping methods, Sanger DNA sequencing, pyrosequencing, high throughput sequencing, targeted sequencing using circularizing probes, targeted sequencing using capture by hybridization probes, reversible dye terminator sequencing, sequencing by ligation, sequencing by hybridization, other methods of DNA sequencing, other high throughput genotyping platforms, fluorescent in situ hybridization (FISH), comparative genomic hybridization (CGH), array CGH, and multiples or combinations thereof.

In an embodiment, making genotypic measurements may be done on genetic material that is amplified and/or preferentially enriched prior to being measured using a technique or technology that is selected from the group consisting of: Polymerase Chain Reaction (PCR), ligand mediated PCR, degenerative oligonucleotide primer PCR, targeted amplification, PCR, mini-PCR, universal PCR amplification, Multiple Displacement Amplification (MDA), allele-specific PCR, allele-specific amplification techniques, linear amplification methods, ligation of substrate DNA followed by another method of amplification, bridge amplification, padlock probes, circularizing probes, capture by hybridization probes, and combinations thereof.

In an embodiment, the method may also include generating a report comprising the established paternity of the fetus. In an embodiment, the invention may comprise a report disclosing the established paternity of the fetus generated using a method described herein.

Disclosed herein are methods for determining the fraction of DNA originating from a target individual is present in a mixture of DNA that contains DNA from the target individual, and also DNA from at least one other individual. According to aspects illustrated herein, in an embodiment, a method for determining a fraction of DNA from a target individual present in a mixed sample of DNA that comprises DNA from the target individual and DNA from a second individual may include making genotypic measurements at a plurality of polymorphic loci from the mixed sample of DNA, obtaining genotypic data at the plurality of polymorphic loci from the second individual, and determining, on a computer, the fraction of DNA from the target individual present in the mixed sample using the genotypic measurements from the mixed sample of DNA, the genotypic data from the second individual, and probabilistic estimation techniques.

In an embodiment, obtaining genotypic data from the second individual includes making genetic measurements from DNA that consists essentially of DNA from the second individual. In an embodiment, obtaining genotypic data from the second individual may include inferring which genotypic measurements from the genotypic measurements made on the mixed sample of DNA are likely attributable to genetic material from the second individual, and using those genotypic measurements that were inferred to be attributable to genetic material from the second individual as the obtained genotypic measurements.

In an embodiment, inferring the genotypic data of the related individual may also include using allele population frequencies at the loci. In an embodiment, the determined fraction of DNA from a target individual is expressed as a probability of fractions of DNA. In an embodiment, the genotypic measurements made from the mixed sample comprise genotypic measurements made by sequencing the DNA in the mixed sample. In an embodiment, the DNA in the mixed sample is preferentially enriched at the plurality of polymorphic loci prior to making genotypic measurements from the mixed sample of DNA. In an embodiment, the polymorphic loci comprise single nucleotide polymorphisms.

In an embodiment, determining the fraction may also include determining a probability of a plurality of fractions of DNA from the target individual present in the mixed sample of DNA, determining the fraction by selecting the fraction from the plurality of fractions with the largest probability. In an embodiment, determining the fraction may also include determining a probability of a plurality of fractions of DNA from the target individual present in the mixed sample of DNA, using a maximum likelihood estimation technique to determine the most likely fraction, and determining the fraction by selecting the fraction that was determined to be the most likely.

In an embodiment, the target individual is a fetus gestating in a pregnant mother, and the second individual is the pregnant mother. In an embodiment, the method may also include using a platform model that relates genotypic data measured at the polymorphic loci, and using a table that relates maternal genotypes to child genotypes. In an embodiment, the determination also uses genotypic measurements at a plurality of polymorphic loci measured on DNA from the father of the fetus. In an embodiment, the method does not make use of genotypic data from the father of the fetus. In an embodiment, the method does not make use of loci on the Y chromosome. In an embodiment, the invention may comprise a report disclosing an established paternity of the fetus determined using a method disclosed herein for determining the fraction of fetal DNA present in the maternal plasma. In an embodiment, the invention may comprise a report disclosing a ploidy state of the fetus determined using a method disclosed herein for determining the fraction of fetal DNA present in the maternal plasma.

BRIEF DESCRIPTION OF THE DRAWINGS

The presently disclosed embodiments will be further explained with reference to the attached drawings, wherein like structures are referred to by like numerals throughout the several views. The drawings shown are not necessarily to scale, with emphasis instead generally being placed upon illustrating the principles of the presently disclosed embodiments.

FIG. 1 shows the distribution of allele intensities from two parental contexts as measured on maternal plasma.

FIG. 2 shows the distribution of paternity related test statistic for 200 unrelated males and the biological father.

FIG. 3 shows two distributions of intensity ratios for 200 unrelated males and the biological father. Each graph correspond to a different input channels.

FIG. 4 shows the cumulative distribution frequency (cdf) curves for the correlation ratio between the fetal genotypic measurements and the parental genotypic measurements for three cases.

FIG. 5 shows histograms of the correlation ratio between the fetal genotypic measurements and the parental genotypic measurements for three cases.

FIG. 6 shows a histogram of the paternity test statistic for 35 samples as compared to an idealized Gaussian distribution of test statistics for 800 unrelated males.

FIG. 7 shows an example of a report disclosing a paternity exclusion.

FIG. 8 shows an example of a report disclosing a paternity inclusion.

FIG. 9 shows an example of a report disclosing an indeterminate result.

While the above-identified drawings set forth presently disclosed embodiments, other embodiments are also contemplated, as noted in the discussion. This disclosure presents illustrative embodiments by way of representation and not limitation. Numerous other modifications and embodiments can be devised by those skilled in the art which fall within the scope and spirit of the principles of the presently disclosed embodiments.

DETAILED DESCRIPTION

According to aspects illustrated herein, a method is provided for determining whether or not an alleged father is the biological father of a fetus that is gestating in a pregnant mother. In an embodiment, the method includes obtaining genetic material from the alleged father, and obtaining a blood sample from the pregnant mother. In an embodiment, the method may include making genotypic measurements of the alleged father and the pregnant mother, and making genotypic measurements on the free floating DNA (ffDNA, i.e. cfDNA) found in the plasma of the pregnant mother. In an embodiment, the method includes obtaining genotypic data for a set of SNPs of the mother and alleged father of the fetus; making genotypic measurements for the set of SNPs on a mixed sample that comprises DNA from the target individual and also DNA from the mother of the target individual. In an embodiment, the method may include using the genotypic measurements to determine, on a computer, the probability that the alleged father is the biological father of the fetus gestating in the pregnant mother. In an embodiment, the method may include using the genotypic data of the pregnant mother and the alleged father to determine an expected allelic distribution for the genotypic measurements of the fetal/maternal DNA mixture if the alleged father were the biological father of the fetus. In an embodiment, the method may include using the genotypic data of the pregnant mother and genotypic data of a plurality of individuals known not to be the father to determine an expected allelic distribution for the genotypic measurements of the fetal/maternal DNA mixture if the alleged father is not the biological father of the fetus. In an embodiment, the method may involve calculating the probabilities that the alleged father is the biological father of the fetus given the expected allelic distributions, and the actual maternal plasma DNA measurements. In an embodiment, the series of steps outlined in the method results in a transformation of the genetic material of the pregnant mother and the alleged father to produce and determine the correct identity of the biological father of a gestating fetus prenatally and in a non-invasive manner. In an embodiment, determining the likelihood that the alleged father is the biological father includes calling or establishing the alleged father as the biological father if the likelihood that the father is excluded from the allelic distribution created using the plurality of unrelated individuals is above a threshold. In an embodiment, determining the likelihood that the alleged father is the biological father includes calling the alleged father as not the biological father if the likelihood that the alleged father is excluded from the allelic distribution created using the plurality of unrelated individuals is below a threshold. In an embodiment, the paternity determination is made by initially assuming that the alleged father is in fact the father of the child; if the alleged father is incorrect, the child genotypes will not fit these predictions, and the initial assumption is considered to be wrong; however, if the child genotypes do fit the predictions, then the assumption is considered to be correct. Thus, the paternity test considers how well the observed ffDNA fits the child genotypes predicted by the alleged father's genotypes. In an embodiment, an electronic or physical report may be generated stating the paternity determination.

In an embodiment, the paternity determination is made using genetic measurements of free-floating DNA (ffDNA) found in maternal blood, and the genotype information from the mother and alleged father. The general method could be applied to measurements of ffDNA using a variety of platforms such as SNP microarrays, untargeted high throughput sequencing, or targeted sequencing. The methods discussed here address the fact that free-floating fetal DNA is found in maternal plasma at low yet unknown concentrations and is difficult to detect. The paternity test may comprise evaluating the ffDNA measurements and how likely they are to have been generated by the alleged father, based on his genotypes. Regardless of the measurement platform, the test may be based on the genotypes measured at polymorphic locations. In some embodiments the possible alleles at each polymorphic locus may be generalized to A and B, and optionally C, D, and/or E, etc.

In an embodiment, this method involves using allele measurement data from a plurality of loci. In an embodiment, the loci are polymorphic. In an embodiment, some or most of the loci are polymorphic. In an embodiment, the polymorphic loci are single nucleotide polymorphisms (SNPs). In an embodiment, some or most of the polymorphic loci are heterozygous. In an embodiment, it is not necessary to determine which loci are heterozygous in advance of the testing.

In an embodiment, a method disclosed herein uses selective enrichment techniques that preserve the relative allele frequencies that are present in the original sample of DNA at each polymorphic locus from a set of polymorphic loci. In some embodiments the amplification and/or selective enrichment technique may involve PCR techniques such as mini-PCR or ligation mediated PCR, fragment capture by hybridization, or circularizing probes such as Molecular Inversion Probes. In some embodiments, methods for amplification or selective enrichment may involve using PCR primers or other probes where, upon correct hybridization to the target sequence, the 3-prime end or 5-prime end of a nucleotide probe is separated from the polymorphic site of the allele by a small number of nucleotides. In an embodiment, probes in which the hybridizing region is designed to hybridize to a polymorphic site are excluded. These embodiments are improvements over other methods that involve targeted amplification and/or selective enrichment in that they better preserve the original allele frequencies of the sample at each polymorphic locus, whether the sample is pure genomic sample from a single individual or mixture of individuals.

In an embodiment, a method disclosed herein uses highly efficient highly multiplexed targeted PCR to amplify DNA followed by high throughput sequencing to determine the allele frequencies at each target locus. One technique that allows highly multiplexed targeted PCR to perform in a highly efficient manner involves designing primers that are unlikely to hybridize with one another. The PCR probes may be selected by creating a thermodynamic model of potentially adverse interactions, or unintended interactions, between at least 500, at least 1,000, at least 5,000, at least 10,000, at least 20,000, at least 50,000, or at least 100,000 potential primer pairs, or between primers and sample DNA, and then using the model to eliminate designs that are incompatible with other the designs in the pool, or with the sample DNA. Another technique that allows highly multiplexed targeted PCR to perform in a highly efficient manner is using a partial or full nesting approach to the targeted PCR. Using one or a combination of these approaches allows multiplexing of at least 300, at least 800, at least 1,200, at least 4,000 or at least 10,000 primers in a single pool with the resulting amplified DNA comprising a majority of DNA molecules that, when sequenced, will map to targeted loci. Using one or a combination of these approaches allows multiplexing of a large number of primers in a single pool with the resulting amplified DNA comprising greater than 50%, greater than 80%, greater than 90%, greater than 95%, greater than 98%, or greater than 99% DNA molecules that map to targeted loci.

In an embodiment, a method disclosed herein involves determining whether the distribution of observed allele measurements is indicative of a paternity inclusion or exclusion using a maximum likelihood estimation (MLE) technique. The use of a maximum likelihood estimation technique is different from and a significant improvement over methods that use single hypothesis rejection technique in that the resultant determinations will be made with significantly higher accuracy. One reason is that single hypothesis rejection techniques does not contain information on the alternative hypothesis. Another reason is that the maximum likelihood technique allows for the determination of optimal cutoff thresholds for each individual sample. Another reason is that the use of a maximum likelihood technique allows the calculation of a confidence for each paternity determination. The ability to make a confidence calculation for each determination allows a practitioner to know which calls are accurate, and which are more likely to be wrong. In some embodiments, a wide variety of methods may be combined with a maximum likelihood estimation technique to enhance the accuracy of the ploidy calls. In an embodiment, a method disclosed herein involves estimating the fetal fraction of DNA in the mixed sample and using that estimation to calculate both the paternity call (determination) and the confidence of the paternity call.

In an embodiment, the method involves calculating a test statistic that is indicative of the degree of relatedness between a first individual and a second individual, given genotypic measurements at a plurality of polymorphic loci for the first individual, and genotypic measurements at a plurality of polymorphic loci for a mixture of DNA where the mixture of DNA comprises DNA from the second individual and a related individual. In an embodiment, the first individual is an alleged father, the second individual is a gestating fetus, and the related individual is the mother of the fetus. The test statistic may be calculated for the fetus, the mother, and a plurality of individuals known to be unrelated to the fetus, thereby generating a distribution of the metric for unrelated individuals. The test statistic may also be calculated for the fetus, the mother, and the alleged father. A single hypothesis rejection test may be used to determine if the test statistic calculated using the alleged father's genotypic data is part of the distribution of test statistics calculated using the genotypic data of the unrelated individuals. If the test statistic calculated using the alleged father's genotypic data is found to be part of the distribution of test statistics calculated using the genotypic data of the unrelated individuals, then paternity can be excluded, that is, the alleged father may be determined to not be related to the fetus. If the test statistic calculated using the alleged father's genotypic data is found not to be part of the distribution of test statistics calculated using the genotypic data of the unrelated individuals, then paternity can be included, that is, the alleged father may be determined to be related to the fetus.

In an embodiment, the paternity determination involves determining the probability of the measured genotypic data given two possible hypotheses: the hypothesis that the alleged father is the biological father of the fetus, and the hypothesis that the alleged father is not the biological father of the fetus. A probability can then be calculated for each of the hypotheses given the data, and the paternity may be established based on the likelihood of each of the two hypotheses. The determination may utilize genetic measurements made on the maternal plasma, genetic measurements made on DNA from the alleged father, and optionally maternal genotypic data. In an embodiment, the maternal genotypic data can be inferred from the genotypic measurements made on the maternal plasma. In an embodiment, the probability can be determined using a partition of the range of possible fetal fractions; the range of fetal fractions could be anywhere from 0.01% to 99.9%, and the mesh may have increments ranging from 10% to 1%, from 1% to 0.1%, and lower than 0.1%. In an embodiment, the partition of possible fetal fractions may be from 2% to 30%, and the increments are about 1%. In an embodiment, the mesh could be continuous, and the likelihoods could be intergrated over the ranges rather than combined. In an embodiment, the probability can be determined using only one fetal fraction, where that fetal fraction may be determined using any appropriate method. For each possible fetal fraction in the mesh, one can calculate the probability of the data given the two hypotheses. For the hypothesis that the alleged father is the biological father, the alleged father genotypes may be used in the calculation of the probability, while for the hypothesis that the alleged father is not the biological father, population based allele frequency data may additionally be used in the calculation of the probability. In an embodiment, one can use the parent contexts and a platform model in calculating the likelihood of data given hypothesis. In an embodiment, the likelihoods can be combined over all fetal fractions in the partition, over all mother genotypes, and over all father genotypes. In an embodiment, the parental genotypes may be probabilistic (e.g. at a given SNP, a parent may have the genotype GT with 99% chance, GG with 0.5% chance, and TT with 0.5% chance; in another embodiment the parental genotypes may take on one value (e.g. at a given SNP, a parent has the genotype GT). In some embodiments the terms probability and likelihood may be interchangeable, as in common parlance; in other embodiments, the two terms may not be interchangeable, and may be read as one skilled in the art in statistics would read them.

In some methods known in the art, fetal fraction is determined using measurements made at loci that are found exclusively on the paternal genotype, for example, loci that are found exclusively on the Y-chromosome, or the Rhesus-D gene. Unfortunately, these methods require either that the fetus is male (in the case where the loci are found exclusively on the Y chromosome) or that a gene or set of genes can be identified prior to measurements of the DNA where those genes are present on the paternal genotype, and not present in the maternal genotype. An additional complication is that in the context of paternity testing, it is not known whether or not the alleged father is the biological father, and therefore, with the exception of the Y-chromosome specific loci, it is not possible to determine what loci may be present on the father, and not on the mother. Therefore, in the context of paternity testing, it is not currently possible to determine the fetal fraction when the fetus is a female, and when the fetus is male, fetal fraction can only be determined by using Y-chromosome specific loci. In an embodiment, a method is disclosed herein for determining the fraction of fetal DNA that is present in the mixture of DNA comprising maternal and fetal DNA. In an embodiment, the method can determine the fetal fraction of fetal DNA that is present in the mixture of DNA comprising maternal and fetal DNA using genotypic measurements from autosomal chromosomes. In an embodiment, the method can determine the fetal fraction of fetal DNA that is present in the mixture of DNA comprising maternal and fetal DNA irrespective of the sex of the fetus. In an embodiment, the method can determine the fetal fraction of fetal DNA that is present in the mixture of DNA comprising maternal and fetal DNA irrespective of what genes the mother and alleged father may have. The instant method does not require that the fetus be male, or that a locus or loci can be identified that are present on the father and not on the mother. The instant method does not require that the paternal genotype be known. The instant method does not require that the maternal genotype be known, as it can be inferred from the measurements made on the DNA in the maternal plasma, which comprises a mixture of both fetal and maternal DNA.

In an embodiment, the distribution of polymorphic loci can be modeled using a binomial distribution. In an embodiment, the distribution of polymorphic loci can be modeled using a beta-binomial distribution. By using the beta-binomial distribution as a model for allele distribution, one can more accurately model likely allele measurements than when using other distributions; this can result in more accurate paternity determinations.

In an embodiment, a method disclosed herein takes into account the tendency for the data to be noisy and contain errors by attaching a probability to each measurement. The use of maximum likelihood techniques to choose the correct hypothesis from the set of hypotheses that were made using the measurement data with attached probabilistic estimates makes it more likely that the incorrect measurements will be discounted, and the correct measurements will be used in the calculations that lead to the paternity determination. To be more precise, this method systematically reduces the influence of incorrectly measured data on the paternity determination. This is an improvement over methods where all data is assumed to be equally correct or methods where outlying data is arbitrarily excluded from calculations leading to a paternity determination. In an embodiment, individual SNPs are weighted by expected measurement variance based on the SNP quality and observed depth of read; this may result in an increase in the accuracy of the resulting statistic, resulting in an increase of the accuracy of the paternity call significantly, especially in borderline cases.

The methods described herein are particularly advantageous when used on samples where a small amount of DNA is available, or where the percent of fetal DNA is low. This is due to the correspondingly higher allele dropout rate that may occur when only a small amount of DNA is available and/or the correspondingly higher fetal allele dropout rate when the percent of fetal DNA is low in a mixed sample of fetal and maternal DNA. A high allele dropout rate, meaning that a large percentage of the alleles were not measured for the target individual, results in poorly accurate fetal fraction calculations, and poorly accurate paternity determinations. The methods described herein allow for an accurate ploidy determination to be made when the percent of molecules of DNA that are fetal in the mixture is less than 40%, less than 30%, less than 20%, less than 10%, less than 8%, less than 6%, less than 4%, and even less than 3%.

In an embodiment, it is possible to determine the paternity of an individual based on measurements when that individual's DNA is mixed with DNA of a related individual. In an embodiment, the mixture of DNA is the free floating DNA found in maternal plasma, which may include DNA from the mother, with known genotype, and which may be mixed with DNA of the fetus, with unknown genotype. The paternity of the fetus can then be determined by looking at the actual measurements, and determining the likelihood of paternity given the observed data. In some embodiments, a method disclosed herein could be used in situations where there is a very small amount of DNA present, such as in forensic situations, where one or a few cells are available (typically less than ten cells, less than twenty cells, less than 40 cells, less than 100 cells, or an equivalent amount of DNA.) In some embodiments, a method disclosed herein could be used in situations where the DNA is highly fragmented, such as ffDNA found in plasma. In these embodiments, a method disclosed herein serves to make paternity calls from a small amount of DNA that is not contaminated by other DNA, but where the paternity calling very difficult due to the small amount of DNA. The genetic measurements used as part of these methods could be made on any sample comprising DNA or RNA, for example but not limited to: blood, plasma, body fluids, urine, hair, tears, saliva, tissue, skin, fingernails, blastomeres, embryos, amniotic fluid, chorionic villus samples, feces, bile, lymph, cervical mucus, semen, or other cells or materials comprising nucleic acids. In an embodiment, a method disclosed herein could be run with nucleic acid detection methods such as sequencing, microarrays, qPCR, digital PCR, or other methods used to measure nucleic acids. In some embodiments, a method disclosed herein involves calculating, on a computer, allele ratios at the plurality of polymorphic loci from the DNA measurements made on the processed samples. In some embodiments, a method disclosed herein involves calculating, on a computer, allele ratios or allelic distributions at a plurality of polymorphic loci from the DNA measurements made on the processed samples along with any combination of other improvements described in this disclosure.

Further discussion of these points may be found elsewhere in this document.

Non-Invasive Prenatal Paternity Testing (NPPT)

The process of non-invasive prenatal paternity testing involves a number of steps. Some of the steps may include: (1) obtaining the genetic material from the fetus; (2) enriching the genetic material of the fetus that may be in a mixed sample, ex vivo; (3) amplifying the genetic material, ex vivo; (4) preferentially enriching specific loci in the genetic material, ex vivo; (5) measuring the genetic material, ex vivo; and (6) analyzing the genotypic data, on a computer, and ex vivo. Methods to reduce to practice these six and other relevant steps are described herein. At least some of the method steps are not directly applied on the body. In an embodiment, the present disclosure relates to methods of treatment and diagnosis applied to tissue and other biological materials isolated and separated from the body. At least some of the method steps are executed on a computer.

Some embodiments of the present disclosure allow a clinician to determine the genetic state of a fetus, specifically its biological relationship to another individual, that is gestating in a mother in a non-invasive manner such that the health of the baby is not put at risk by the collection of the genetic material of the fetus, and that the mother is not required to undergo an invasive procedure.

Modem technological advances have resulted in the ability to measure large amounts of genetic information from a genetic sample using such methods as high throughput sequencing and genotyping arrays. The methods disclosed herein allow a clinician to take greater advantage of the large amounts of data available, and make a more accurate diagnosis of the fetal genetic identity. In an embodiment, an informatics based method may result in paternity determinations of higher accuracy than by methods currently known in the art. The details of a number of embodiments are given below. Different embodiments may involve different combinations of the aforementioned steps. Various combinations of the different embodiments of the different steps may be used interchangeably.

In an embodiment, a blood sample is taken from a pregnant mother, and the free floating DNA in the plasma of the mother's blood, which contains a mixture of both DNA of maternal origin, and DNA of fetal origin, is isolated and used to determine the ploidy status of the fetus. In an embodiment, a method disclosed herein involves preferential enrichment of those DNA sequences in a mixture of DNA that correspond to polymorphic alleles in a way that the allele ratios and/or allele distributions remain reasonably consistent upon enrichment. In an embodiment, the method involves amplifying the isolated DNA using whole genome amplification (WGA). In an embodiment, a method disclosed herein involves targeted PCR based amplification such that a high percentage of the resulting molecules correspond to targeted loci. In an embodiment, a method disclosed herein involves sequencing a mixture of DNA that contains both DNA of maternal origin, and DNA of fetal origin. In an embodiment, the method involves measuring the amplified DNA using a microarray designed to detect nucleic acid sequences such as a SNP array. In an embodiment, a method disclosed herein involves using measured allele distributions to determine the paternity of a fetus that is gestating in a mother. In an embodiment, a method disclosed herein involves reporting the determined paternity state to a clinician. In an embodiment, a method disclosed herein involves taking a clinical action, for example, performing follow up invasive testing such as chorionic villus sampling or amniocentesis, preparing for the birth of a child, or an elective termination of a fetus.

This application makes reference to U.S. Utility application Ser. No. 11/603,406, filed Nov. 28, 2006 (US Publication No.: 20070184467); U.S. Utility application Ser. No. 12/076,348, filed Mar. 17, 2008 (US Publication No.: 20080243398); PCT Utility Application Serial No. PCT/US09/52730, filed Aug. 4, 2009 (PCT Publication No.: WO/2010/017214); PCT Utility Application Serial No. PCT/US10/050824, filed Sep. 30, 2010 (PCT Publication No.: WO/2011/041485), U.S. Utility application Ser. No. 13/110,685, filed May 18, 2011, and U.S. Utility application Ser. No. 13/300,235, filed Nov. 18, 2011. Some of the vocabulary used in this filing may have its antecedents in these references. Some of the concepts described herein may be better understood in light of the concepts found in these references.

Screening Maternal Blood Comprising Free Floating Fetal DNA

The methods described herein may be used to help determine whether a child, fetus, or other target individual is genetically related to another individual. In some embodiment, this may be done in cases where the genetic material of the target individual is found in the presence of a quantity of genetic material from another individual. In one embodiment, the method may be used to help determine whether a fetus is genetically related to an alleged father using the free floating fetal DNA found in the maternal blood, along with a genetic sample from the father and optionally the mother. In an embodiment, the fetus may have originated from an egg from an egg donor such that the fetus is not genetically related to the mother in which the fetus is gestating. In an embodiment, the method may be applicable in cases where the amount of target DNA is in any proportion with the non-target DNA; for example, the target DNA could make up anywhere between 0.000001 and 99.999999% of the DNA present. In an embodiment, the non-target contaminating DNA could be from a plurality of individuals; it is advantageous where genetic data from some or all of the relevant non-target individual(s) is known, or where genetic samples from said related individuals are available. In an embodiment, a method disclosed herein can be used to determine genotypic data of a fetus from maternal blood that contains fetal DNA. It may also be used in a case where there are multiple fetuses in the uterus of a pregnant woman, or where other contaminating DNA may be present in the sample, for example from other already born siblings.

This technique may make use of the phenomenon of fetal blood cells gaining access to maternal circulation through the placental villi. Ordinarily, only a very small number of fetal cells enter the maternal circulation in this fashion (not enough to produce a positive Kleihauer-Betke test for fetal-maternal hemorrhage). The fetal cells can be sorted out and analyzed by a variety of techniques to look for particular DNA sequences, but without the risks that invasive procedures inherently have. This technique may also make use of the phenomenon of free floating fetal DNA gaining access to maternal circulation by DNA release following apoptosis of placental tissue where the placental tissue in question contains DNA of the same genotype as the fetus. The free floating DNA found in maternal plasma has been shown to contain fetal DNA in proportions as high as 30-40% fetal DNA.

In an embodiment, blood may be drawn from a pregnant woman. Research has shown that maternal blood may contain a small amount of free floating DNA from the fetus, in addition to free floating DNA of maternal origin. In addition, there also may be nucleated fetal blood cells comprising DNA of fetal origin, in addition to many blood cells of maternal origin, which typically do not contain nuclear DNA. There are many methods know in the art to isolate fetal DNA, or create fractions enriched in fetal DNA. For example, chromatography has been show to create certain fractions that are enriched in fetal DNA.

Once the sample of maternal blood, plasma, or other fluid, drawn in a relatively non-invasive manner, and that contains an amount of fetal DNA, either cellular or free floating, either enriched in its proportion to the maternal DNA, or in its original ratio, is in hand, one may genotype the DNA found in said sample. In some embodiments, the blood may be drawn using a needle to withdraw blood from a vein, for example, the basilica vein. The method described herein can be used to determine genotypic data of the fetus. For example, it can be used to determine the ploidy state at one or more chromosomes, it can be used to determine the identity of one or a set of SNPs, including insertions, deletions, and translocations. It can be used to determine one or more haplotypes, including the parent of origin of one or more genotypic features. It can also be used to determine the degree of relatedness between the fetus an another individual.

Note that this method will work with any nucleic acids that can be used for any genotyping and/or sequencing methods, such as the ILLUMINA INFINIUM ARRAY platform, AFFYMETRIX GENECHIP, ILLUMINA GENOME ANALYZER, or LIFE TECHNOLOGIES' SOLID SYSTEM, along with the genotypic data measured therefrom. This includes extracted free-floating DNA from plasma or amplifications (e.g. whole genome amplification, PCR) of the same; genomic DNA from other cell types (e.g. human lymphocytes from whole blood) or amplifications of the same. For preparation of the DNA, any extraction or purification method that generates genomic DNA suitable for the one of these platforms will work as well. This method could work equally well with samples of RNA. In an embodiment, storage of the samples may be done in a way that will minimize degradation (e.g. below freezing, at about −20 C, or at a lower temperature).

Parental Support

Some embodiments may be used in combination with the PARENTAL SUPPORT™ (PS) method, embodiments of which are described in U.S. application Ser. No. 11/603,406 (US Publication No.: 20070184467), U.S. application Ser. No. 12/076,348 (US Publication No.: 20080243398), U.S. application Ser. No. 13/110,685, PCT Application PCT/US09/52730 (PCT Publication No.: WO/2010/017214), PCT Application No. PCT/US10/050824 (PCT Publication No.: WO/2011/041485), PCT Application No. PCT/US2011/037018 (PCT Publication No.: WO/2011/146632), and PCT Application No. PCT/US2011/61506, which are incorporated herein by reference in their entirety. PARENTAL SUPPORT™ is an informatics based approach that can be used to analyze genetic data. In some embodiments, the methods disclosed herein may be considered as part of the PARENTAL SUPPORT™ method. In some embodiments, The PARENTAL SUPPORT™ method is a collection of methods that may be used to determine the genetic data of a target individual, with high accuracy, of one or a small number of cells from that individual, or of a mixture of DNA consisting of DNA from the target individual and DNA from one or a plurality of other individuals, specifically to determine disease-related alleles, other alleles of interest, the ploidy state of one or a plurality of chromosomes in the target individual, and or the extent of relationship of another individual to the target individual. PARENTAL SUPPORT™ may refer to any of these methods. PARENTAL SUPPORT™ is an example of an informatics based method.

The PARENTAL SUPPORT™ method makes use of known parental genetic data, i.e. haplotypic and/or diploid genetic data of the mother and/or the father, together with the knowledge of the mechanism of meiosis and the imperfect measurement of the target DNA, and possibly of one or more related individuals, along with population based crossover frequencies, in order to reconstruct, in silico, the genotype at a plurality of alleles, and/or the paternity state of an embryo or of any target cell(s), and the target DNA at the location of key loci with a high degree of confidence. The PARENTAL SUPPORT™ method makes use of known parental genetic data, i.e. haplotypic and/or diploid genetic data of the mother and/or the father, together with the knowledge of the mechanism of meiosis and the imperfect measurement of the target DNA, to create hypotheses about what genetic data may be expected for different situations, to calculate the likelihood of each of the situations given the observed genetic data, thereby determining which situation is most likely. In some embodiments the situation is question may include whether the target individual has inherited a disease linked haplotype of interest, whether the target individual has inherited a phenotype linked haplotype of interest, whether the target individual has one or more aneuploid chromosomes, and/or whether the target individual is related to an individual of interest, and what the degree of relationship may be. The PARENTAL SUPPORT™ method allows the cleaning of noisy genetic data. PARENTAL SUPPORT™ may be particularly relevant where only a small fraction of the genetic material available is from the target individual (e.g. NPD or NPPT) and where direct measurements of the genotypes are inherently noisy due to the contaminating DNA signal from another individual. The PARENTAL SUPPORT™ method is able to reconstruct highly accurate ordered diploid allele sequences on the embryo, together with copy number of chromosomes segments, even though the conventional, unordered diploid measurements may be characterized by high rates of allele dropouts, drop-ins, variable amplification biases and other errors. The method may employ both an underlying genetic model and an underlying model of measurement error. The genetic model may determine both allele probabilities at each SNP and crossover probabilities between SNPs. Allele probabilities may be modeled at each SNP based on data obtained from the parents and model crossover probabilities between SNPs based on data obtained from the HapMap database, as developed by the International HapMap Project. Given the proper underlying genetic model and measurement error model, maximum a posteriori (MAP) estimation may be used, with modifications for computationally efficiency, to estimate the correct, ordered allele values at each SNP in the embryo.

Definitions

-   Single Nucleotide Polymorphism (SNP) refers to a single nucleotide     that may differ between the genomes of two members of the same     species. The usage of the term should not imply any limit on the     frequency with which each variant occurs. -   Sequence refers to a DNA sequence or a genetic sequence. It may     refer to the primary, physical structure of the DNA molecule or     strand in an individual. It may refer to the sequence of nucleotides     found in that DNA molecule, or the complementary strand to the DNA     molecule. It may refer to the information contained in the DNA     molecule as its representation in silico. -   Locus refers to a particular region of interest on the DNA of an     individual, which may refer to a SNP, the site of a possible     insertion or deletion, or the site of some other relevant genetic     variation. Disease-linked SNPs may also refer to disease-linked     loci. -   Polymorphic Allele, also “Polymorphic Locus,” refers to an allele or     locus where the genotype varies between individuals within a given     species. Some examples of polymorphic alleles include single     nucleotide polymorphisms, short tandem repeats, deletions,     duplications, and inversions. -   Polymorphic Site refers to the specific nucleotides found in a     polymorphic region that vary between individuals. -   Allele refers to the genes that occupy a particular locus. -   Genetic Data also “Genotypic Data” refers to the data describing     aspects of the genome of one or more individuals. It may refer to     one or a set of loci, partial or entire sequences, partial or entire     chromosomes, or the entire genome. It may refer to the identity of     one or a plurality of nucleotides; it may refer to a set of     sequential nucleotides, or nucleotides from different locations in     the genome, or a combination thereof. Genotypic data is typically in     silico, however, it is also possible to consider physical     nucleotides in a sequence as chemically encoded genetic data.     Genotypic Data may be said to be “on,” “of,” “at,” “from” or “on”     the individual(s). Genotypic Data may refer to output measurements     from a genotyping platform where those measurements are made on     genetic material. -   Genetic Material also “Genetic Sample” refers to physical matter,     such as tissue or blood, from one or more individuals comprising DNA     or RNA -   Confidence refers to the statistical likelihood that the called SNP,     allele, set of alleles, ploidy call, or paternity call is correct. -   Aneuploidy refers to the state where the wrong number of chromosomes     is present in a cell. In the case of a somatic human cell it may     refer to the case where a cell does not contain 22 pairs of     autosomal chromosomes and one pair of sex chromosomes. In the case     of a human gamete, it may refer to the case where a cell does not     contain one of each of the 23 chromosomes. In the case of a single     chromosome type, it may refer to the case where more or less than     two homologous but non-identical chromosome copies are present, or     where there are two chromosome copies present that originate from     the same parent. -   Chromosome may refer to a single chromosome copy, meaning a single     molecule of DNA of which there are 46 in a normal somatic cell; an     example is ‘the maternally derived chromosome 18’. Chromosome may     also refer to a chromosome type, of which there are 23 in a normal     human somatic cell; an example is ‘chromosome 18’. -   Monosomy refers to the state where a cell only contains one of a     chromosome type. -   Disomy refers to the state where a cell contains two of a chromosome     type. -   Uniparental Disomy refers to the state where a cell contains two of     a chromosome type, and where both chromosomes originate from one     parent. -   Trisomy refers to the state where a cell contains three of a     chromosome type. -   The State of the Genetic Material or simply “Genetic State” may     refer to the identity of a set of SNPs on the DNA, to the phased     haplotypes of the genetic material, or to the sequence of the DNA,     including insertions, deletions, repeats and mutations. It may also     refer to the ploidy state of one or more chromosomes, chromosomal     segments, or set of chromosomal segments. -   Establishing the Paternity or “Determining the Paternity” refers to     establishing or determining that an alleged father either is or is     not the biological father of a gestating fetus, or determining or     establishing the likelihood that an alleged father is the biological     father of the fetus. -   Paternity Determination refers to the determination that the alleged     father is or is not the biological father of the fetus. A paternity     determination is the result of establishing, calling or determining     the paternity. -   Paternity refers to the identity of the biological father of an     individual. -   Paternity Inclusion refers to establishing that an alleged father is     the biological father of a fetus. -   Paternity Exclusion refers to establishing that an alleged father is     not the biological father of a fetus. -   Alleged Father refers to a male whose paternal relationship to a     fetus is in question. -   Biological Father of an individual refers to the male whose genetic     material was inherited by the individual. -   Allelic Ratio refers to the ratio between the amount of each allele     at a polymorphic locus that is present in a sample or in an     individual. When the sample is measured by sequencing, the allelic     ratio may refer to the ratio of sequence reads that map to each     allele at the locus. When the sample is measured by an intensity     based measurement method, the allele ratio may refer to the ratio of     the amounts of each allele present at that locus as estimated by the     measurement method. -   Allelic Distribution, or ‘allele count distribution’ refers to the     relative amount of each allele that is present for each locus in a     set of loci. An allelic distribution can refer to an individual, to     a sample, or to a set of measurements made on a sample. In the     context of sequencing, the allelic distribution refers to the number     or probable number of reads that map to a particular allele for each     allele in a set of polymorphic loci. The allele measurements may be     treated probabilistically, that is, the likelihood that a given     allele is present for a give sequence read is a fraction between 0     and 1, or they may be treated in a binary fashion, that is, any     given read is considered to be exactly zero or one copies of a     particular allele. -   Allelic Bias refers to the degree to which the measured ratio of     alleles at a heterozygous locus is different to the ratio that was     present in the original sample of DNA. The degree of allelic bias at     a particular locus is equal to the observed allelic ratio at that     locus, as measured, divided by the ratio of alleles in the original     DNA sample at that locus. Allelic bias may be defined to be greater     than one, such that if the calculation of the degree of allelic bias     returns a value, x, that is less than 1, then the degree of allelic     bias may be restated as 1/x. Allelic bias maybe due to amplification     bias, purification bias, or some other phenomenon that affects     different alleles differently. -   Primer, also “PCR probe” refers to a single DNA molecule (a DNA     oligomer) or a collection of DNA molecules (DNA oligomers) where the     DNA molecules are identical, or nearly so, and where the primer     contains a region that is designed to hybridize to a targeted     polymorphic locus, and may contain a priming sequence designed to     allow PCR amplification. A primer may also contain a molecular     barcode. A primer may contain a random region that differs for each     individual molecule. -   Hybrid Capture Probe refers to any nucleic acid sequence, possibly     modified, that is generated by various methods such as PCR or direct     synthesis and intended to be complementary to one strand of a     specific target DNA sequence in a sample. The exogenous hybrid     capture probes may be added to a prepared sample and hybridized     through a denature-reannealing process to form duplexes of     exogenous-endogenous fragments. These duplexes may then be     physically separated from the sample by various means. -   Sequence Read refers to data representing a sequence of nucleotide     bases that were measured using a clonal sequencing method. Clonal     sequencing may produce sequence data representing single, or clones,     or clusters of one original DNA molecule. A sequence read may also     have associated quality score at each base position of the sequence     indicating the probability that nucleotide has been called     correctly. -   Mapping a sequence read is the process of determining a sequence     read's location of origin in the genome sequence of a particular     organism. The location of origin of sequence reads is based on     similarity of nucleotide sequence of the read and the genome     sequence. -   Homozygous refers to having similar alleles at corresponding     chromosomal loci. -   Heterozygous refers to having dissimilar alleles at corresponding     chromosomal loci. -   Heterozygosity Rate refers to the rate of individuals in the     population having heterozygous alleles at a given locus. The     heterozygosity rate may also refer to the expected or measured ratio     of alleles, at a given locus in an individual, or a sample of DNA. -   Haplotype refers to a combination of alleles at multiple loci that     are typically inherited together on the same chromosome. Haplotype     may refer to as few as two loci or to an entire chromosome depending     on the number of recombination events that have occurred between a     given set of loci. Haplotype can also refer to a set of single     nucleotide polymorphisms (SNPs) on a single chromatid that are     statistically associated. -   Haplotypic Data, also “Phased Data” or “Ordered Genetic Data,”     refers to data from a single chromosome in a diploid or polyploid     genome, i.e., either the segregated maternal or paternal copy of a     chromosome in a diploid genome. -   Phasing refers to the act of determining the haplotypic genetic data     of an individual given unordered, diploid (or polyploid) genetic     data. It may refer to the act of determining which of two genes at     an allele, for a set of alleles found on one chromosome, are     associated with each of the two homologous chromosomes in an     individual. -   Phased Data refers to genetic data where one or more haplotypes have     been determined. -   Fetal refers to “of the fetus,” or “of the region of the placenta     that is genetically similar to the fetus”. In a pregnant woman, some     portion of the placenta is genetically similar to the fetus, and the     free floating fetal DNA found in maternal blood may have originated     from the portion of the placenta with a genotype that matches the     fetus. -   DNA of Fetal Origin refers to DNA that was originally part of a cell     whose genotype was essentially equivalent to that of the fetus. Note     that the genetic information in half of the chromosomes in a fetus     is inherited from the mother of the fetus; in some embodiments, the     DNA from these maternally inherited chromosomes that came from a     fetal cell is considered to be “of fetal origin,” and not “of     maternal origin.” -   DNA of Maternal Origin refers to DNA that was originally part of a     cell whose genotype was essentially equivalent to that of the     mother. -   Child may refer to an embryo, a blastomere, or a fetus. Note that in     the presently disclosed embodiments, the concepts described apply     equally well to individuals who are a born child, a fetus, an embryo     or a set of cells therefrom. The use of the term child may simply be     meant to connote that the individual referred to as the child is the     genetic offspring of the parents. -   Parent refers to the genetic mother or father of an individual. An     individual typically has two parents, a mother and a father, though     this may not necessarily be the case such as in genetic or     chromosomal chimerism. -   Mother may refer to the biological mother of an individual, and/or     it may refer to the women who is carrying the individual as he/she     gestates. -   Parental Context refers to the genetic state of a given SNP, on each     of the two relevant chromosomes for one or both of the two parents     of the target. -   Maternal Plasma refers to the plasma portion of the blood from a     female who is pregnant. -   Clinical Decision refers to any decision to take or not take an     action that has an outcome that affects the health or survival of an     individual. In the context of prenatal paternity testing, a clinical     decision may refer to a decision to abort or not abort a fetus. A     clinical decision may also refer to a decision to conduct further     testing, or to take actions to prepare for the birth of a child. -   Diagnostic Box refers to one or a combination of machines designed     to perform one or a plurality of aspects of the methods disclosed     herein. In an embodiment, the diagnostic box may be placed at a     point of patient care. In an embodiment, the diagnostic box may     perform targeted amplification followed by sequencing. In an     embodiment the diagnostic box may function alone or with the help of     a technician. -   Informatics Based Method or ‘informatics based approach’ refers to a     method that relies heavily on statistics to make sense of a large     amount of data. In the context of prenatal diagnosis, it refers to a     method designed to determine the ploidy state at one or more     chromosomes or the allelic state at one or more alleles by     statistically inferring the most likely state, rather than by     directly physically measuring the state, given a large amount of     genetic data, for example from a molecular array or sequencing. In     an embodiment of the present disclosure, the informatics based     technique may be one disclosed in this patent. In an embodiment of     the present disclosure it may be PARENTAL SUPPORT™. -   Preferential Enrichment of DNA that corresponds to one or a     plurality of loci, or preferential enrichment of DNA at one or a     plurality of loci, refers to any method that results in the     percentage of molecules of DNA in a post-enrichment DNA mixture that     correspond to the loci being higher than the percentage of molecules     of DNA in the pre-enrichment DNA mixture that correspond to the     loci. The method may involve selective amplification of DNA     molecules that correspond to the loci. The method may involve     removing DNA molecules that do not correspond to the loci. -   Amplification refers to a method that increases the number of copies     of a molecule of DNA. -   Selective Amplification may refer to a method that increases the     number of copies of a particular molecule of DNA, or molecules of     DNA that correspond to a particular region of DNA. It may also refer     to a method that increases the number of copies of a particular     targeted molecule of DNA, or targeted region of DNA more than it     increases non-targeted molecules or regions of DNA. Selective     amplification may be a method of preferential enrichment. -   Universal Priming Sequence refers to a DNA sequence that may be     appended to a population of target DNA molecules, for example by     ligation, PCR, or ligation mediated PCR. Once added to the     population of target molecules, primers specific to the universal     priming sequences can be used to amplify the target population using     a single pair of amplification primers. Universal priming sequences     are typically not related to the target sequences. -   Universal Adapters, or ‘ligation adaptors’ or ‘library tags’ are DNA     molecules containing a universal priming sequence that can be     covalently linked to the 5-prime and 3-prime end of a population of     target double stranded DNA molecules. The addition of the adapters     provides universal priming sequences to the 5-prime and 3-prime end     of the target population from which PCR amplification can take     place, amplifying all molecules from the target population, using a     single pair of amplification primers. -   Targeting refers to a method used to selectively amplify or     otherwise preferentially enrich those molecules of DNA that     correspond to a set of loci, in a mixture of DNA. -   Hypothesis refers to the possibility that the alleged father is the     biological father of the fetus, or that the alleged father is not     the biological father of the fetus. -   Determining, establishing, and calculating may be used     interchangeably.

Parental Contexts

The parental context refers to the genetic state of a given allele, on each of the two relevant chromosomes for one or both of the two parents of the target. Note that in an embodiment, the parental context does not refer to the allelic state of the target, rather, it refers to the allelic state of the parents. The parental context for a given SNP may consist of four base pairs, two paternal and two maternal; they may be the same or different from one another. It is typically written as “m₁m₂|f₁f₂” where m₁ and m₂ are the genetic state of the given SNP on the two maternal chromosomes, and f₁ and f₂ are the genetic state of the given SNP on the two paternal chromosomes. In some embodiments, the parental context may be written as “f₁f₂|m₁m₂.” Note that subscripts “1” and “2” refer to the genotype, at the given allele, of the first and second chromosome; also note that the choice of which chromosome is labeled “1” and which is labeled “2” may be arbitrary.

Note that in this disclosure, A and B are often used to generically represent base pair identities; A or B could equally well represent C (cytosine), G (guanine), A (adenine) or T (thymine). For example, if, at a given SNP based allele, the mother's genotype was T at that SNP on one chromosome, and G at that SNP on the homologous chromosome, and the father's genotype at that allele is G at that SNP on both of the homologous chromosomes, one may say that the target individual's allele has the parental context of AB|BB; it could also be said that the allele has the parental context of AB|AA. Note that, in theory, any of the four possible nucleotides could occur at a given allele, and thus it is possible, for example, for the mother to have a genotype of AT, and the father to have a genotype of GC at a given allele. However, empirical data indicate that in most cases only two of the four possible base pairs are observed at a given allele. It is possible, for example when using single tandem repeats, to have more than two parental, more than four and even more than ten contexts. In this disclosure the discussion assumes that only two possible base pairs will be observed at a given allele, although the embodiments disclosed herein could be modified to take into account the cases where this assumption does not hold.

A “parental context” may refer to a set or subset of target SNPs that have the same parental context. For example, if one were to measure 1000 alleles on a given chromosome on a target individual, then the context AA|BB could refer to the set of all alleles in the group of 1,000 alleles where the genotype of the mother of the target was homozygous, and the genotype of the father of the target is homozygous, but where the maternal genotype and the paternal genotype are dissimilar at that locus. If the parental data is not phased, and thus AB=BA, then there are nine possible parental contexts: AA|AA, AA|AB, AA|BB, AB|AA, AB|AB, AB|BB, BB|AA, BB|AB, and BB|BB. If the parental data is phased, and thus AB≠BA, then there are sixteen different possible parental contexts: AA|AA, AA|AB, AA|BA, AA|BB, AB|AA, AB|AB, AB|BA, AB|BB, BA|AA, BA|AB, BA|BA, BA|BB, BB|AA, BB|AB, BB|BA, and BB|BB. Every SNP allele on a chromosome, excluding some SNPs on the sex chromosomes, has one of these parental contexts. The set of SNPs wherein the parental context for one parent is heterozygous may be referred to as the heterozygous context.

Different Implementations of the Presently Disclosed Embodiments

Method are disclosed herein for determining the paternity of a target individual. The target individual may be a blastomere, an embryo, or a fetus. In some embodiments of the present disclosure, a method for determining the paternity of an individual may include any of the steps described in this document, and combinations thereof:

In some embodiments the source of the genetic material to be used in determining the paternity of the fetus may be fetal cells, such as nucleated fetal red blood cells, isolated from the maternal blood. The method may involve obtaining a blood sample from the pregnant mother. In some embodiments of the present disclosure, the genetic material to be used in determining the paternity of the fetus may free floating DNA from maternal plasma, where the free floating DNA may be comprised of a mixture of fetal and maternal DNA.

In some embodiments, the source of the genetic material of the fetus may be fetal cells, such as nucleated fetal red blood cells, isolated from the maternal blood. The method may involve obtaining a blood sample from the pregnant mother. The method may involve isolating a fetal red blood cell using visual techniques, based on the idea that a certain combination of colors are uniquely associated with nucleated red blood cell, and a similar combination of colors is not associated with any other present cell in the maternal blood. The combination of colors associated with the nucleated red blood cells may include the red color of the hemoglobin around the nucleus, which color may be made more distinct by staining, and the color of the nuclear material which can be stained, for example, blue. By isolating the cells from maternal blood and spreading them over a slide, and then identifying those points at which one sees both red (from the Hemoglobin) and blue (from the nuclear material) one may be able to identify the location of nucleated red blood cells. One may then extract those nucleated red blood cells using a micromanipulator, use genotyping and/or sequencing techniques to measure aspects of the genotype of the genetic material in those cells.

In one embodiment, one may stain the nucleated red blood cell with a die that only fluoresces in the presence of fetal hemoglobin and not maternal hemoglobin, and so remove the ambiguity between whether a nucleated red blood cell is derived from the mother or the fetus. Some embodiments of the present disclosure may involve staining or otherwise marking nuclear material. Some embodiments of the present disclosure may involve specifically marking fetal nuclear material using fetal cell specific antibodies.

There are many other ways to isolate fetal cells from maternal blood, or fetal DNA from maternal blood, or to enrich samples of fetal genetic material in the presence of maternal genetic material. Some of these methods are listed here, but this is not intended to be an exhaustive list. Some appropriate techniques are listed here for convenience: using fluorescently or otherwise tagged antibodies, size exclusion chromatography, magnetically or otherwise labeled affinity tags, epigenetic differences, such as differential methylation between the maternal and fetal cells at specific alleles, density gradient centrifugation succeeded by CD45/14 depletion and CD71-positive selection from CD45/14 negative-cells, single or double Percoll gradients with different osmolalities, or galactose specific lectin method.

In some embodiments, the genetic sample may be prepared, isolated and/or purified. In some embodiments, the sample may be centrifuged to separate various layers. In some embodiments the preparation of the DNA may involve amplification, separation, purification by chromatography, purification by electrophoresis, filtration, liquid separation, isolation, precipitation, preferential enrichment, preferential amplification, targeted amplification, or any of a number of other techniques either known in the art or described herein.

In some embodiments, the method of the present disclosure may involve amplifying DNA. Amplification of the DNA, a process which transforms a small amount of genetic material to a larger amount of genetic material that comprises a similar set of genetic data, can be done by a wide variety of methods, including, but not limited to polymerase chain reaction (PCR). One method of amplifying DNA is whole genome amplification (WGA). There are a number of methods available for WGA: ligation-mediated PCR (LM-PCR), degenerate oligonucleotide primer PCR (DOP-PCR), and multiple displacement amplification (MDA). In LM-PCR, short DNA sequences called adapters are ligated to blunt ends of DNA. These adapters contain universal amplification sequences, which are used to amplify the DNA by PCR. In DOP-PCR, random primers that also contain universal amplification sequences are used in a first round of annealing and PCR. Then, a second round of PCR is used to amplify the sequences further with the universal primer sequences. MDA uses the phi-29 polymerase, which is a highly processive and non-specific enzyme that replicates DNA and has been used for single-cell analysis. Single-cell whole genome amplification has been used successfully for a variety of applications for a number of years. There are other methods of amplifying DNA from a sample of DNA. The DNA amplification transforms the initial sample of DNA into a sample of DNA that is similar in the set of sequences, but of much greater quantities. In some cases, amplification may not be required.

In some embodiments, DNA may be amplified using a universal amplification, such as WGA or MDA. In some embodiments, DNA may be amplified by targeted amplification, for example using targeted PCR, or circularizing probes. In some embodiments, the DNA may be preferentially enriched using a targeted amplification method, or a method that results in the full or partial separation of desired from undesired DNA, such as capture by hybridization approaches. In some embodiments, DNA may be amplified by using a combination of a universal amplification method and a preferential enrichment method. A fuller description of some of these methods can be found elsewhere in this document.

The genetic data of the target individual and/or of the related individual can be transformed from a molecular state to an electronic state by measuring the appropriate genetic material using tools and or techniques taken from a group including, but not limited to: genotyping microarrays, and high throughput sequencing. Some high throughput sequencing methods include Sanger DNA sequencing, pyrosequencing, the ILLUMINA SOLEXA platform, ILLUMINA's GENOME ANALYZER, or APPLIED BIOSYSTEM's 454 sequencing platform, HELICOS's TRUE SINGLE MOLECULE SEQUENCING platform, HALCYON MOLECULAR's electron microscope sequencing method, or any other sequencing method. All of these methods physically transform the genetic data stored in a sample of DNA into a set of genetic data that is typically stored in a memory device en route to being processed.

A relevant individual's genetic data may be measured by analyzing substances taken from a group including, but not limited to: the individual's bulk diploid tissue, one or more diploid cells from the individual, one or more haploid cells from the individual, one or more blastomeres from the target individual, extra-cellular genetic material found on the individual, extra-cellular genetic material from the individual found in maternal blood, cells from the individual found in maternal blood, one or more embryos created from a gamete from the related individual, one or more blastomeres taken from such an embryo, extra-cellular genetic material found on the related individual, genetic material known to have originated from the related individual, and combinations thereof.

In some embodiments, the likelihood that an alleged father is the biological father of a fetus may be calculated. In some embodiments, the paternity determination may be used to make a clinical decision. This knowledge, typically stored as a physical arrangement of matter in a memory device, may then be transformed into a report. The report may then be acted upon. For example, the clinical decision may be to terminate the pregnancy; alternately, the clinical decision may be to continue the pregnancy.

In an embodiment of the present disclosure, any of the methods described herein may be modified to allow for multiple targets to come from same target individual, for example, multiple blood draws from the same pregnant mother. This may improve the accuracy of the model, as multiple genetic measurements may provide more data with which the target genotype may be determined. In an embodiment, one set of target genetic data served as the primary data which was reported, and the other served as data to double-check the primary target genetic data. In an embodiment, a plurality of sets of genetic data, each measured from genetic material taken from the target individual, are considered in parallel, and thus both sets of target genetic data serve to help determine the paternity of the fetus.

In an embodiment, the method may be used for the purpose of paternity testing. For example, given the SNP-based genotypic information from the mother, and from a man who may or may not be the genetic father, and the measured genotypic information from the mixed sample, it is possible to determine if the genotypic information of the male indeed represents that actual genetic father of the gestating fetus. A simple way to do this is to simply look at the contexts where the mother is AA, and the possible father is AB or BB. In these cases, one may expect to see the father contribution half (AA|AB) or all (AA|BB) of the time, respectively. Taking into account the expected ADO, it is straightforward to determine whether or not the fetal SNPs that are observed are correlated with those of the possible father. Other methods for making a paternity determination are described elsewhere in this document.

In an embodiment of the present disclosure, a pregnant mother would like to determine if a man is the biological father of her fetus. She goes to her doctor, and gives a sample of her blood, and she and her husband gives samples of their own DNA from cheek swabs. A laboratory researcher genotypes the parental DNA using the MDA protocol to amplify the parental DNA, and ILLUMINA INFINIUM arrays to measure the genetic data of the parents at a large number of SNPs. The researcher then spins down the blood, takes the plasma, and isolates a sample of free-floating DNA using size exclusion chromatography. Alternately, the researcher uses one or more fluorescent antibodies, such as one that is specific to fetal hemoglobin to isolate a nucleated fetal red blood cell. The researcher then takes the isolated or enriched fetal genetic material and amplifies it using a library of 70-mer oligonucleotides appropriately designed such that two ends of each oligonucleotide corresponded to the flanking sequences on either side of a target allele. Upon addition of a polymerase, ligase, and the appropriate reagents, the oligonucleotides underwent gap-filling circularization, capturing the desired allele. An exonuclease was added, heat-inactivated, and the products were used directly as a template for PCR amplification. The PCR products were sequenced on an ILLUMINA GENOME ANALYZER. The sequence reads were used as input for the PARENTAL SUPPORT™ method, which then predicted the ploidy state of the fetus. The method determines that the alleged father is not the biological father of the fetus, and calculates a confidence on the determination of 99.98%. A report is generated disclosing both the paternity determination and the confidence of the determination.

In another embodiment a woman who is pregnant wants to know if a man is the biological father of her fetus. The obstetrician takes a blood draw from the mother and father. The blood is sent to a laboratory, where a technician centrifuges the maternal sample to isolate the plasma and the buffy coat. The DNA in the buffy coat and the paternal blood sample are transformed through amplification and the genetic data encoded in the amplified genetic material is further transformed from molecularly stored genetic data into electronically stored genetic data by running the genetic material on a high throughput sequencer to measure the parental genotypes. The plasma sample is preferentially enriched at a set of loci using a 5,000-plex hemi-nested targeted PCR method. The mixture of DNA fragments is prepared into a DNA library suitable for sequencing. The DNA is then sequenced using a high throughput sequencing method, for example, the ILLUMINA GAIIx GENOME ANALYZER. The sequencing transforms the information that is encoded molecularly in the DNA into information that is encoded electronically in computer hardware. An informatics based technique that includes the presently disclosed embodiments, such as PARENTAL SUPPORT™, may be used to determine the paternity of the fetus. This may involve calculating, on a computer, allele counts at the plurality of polymorphic loci from the DNA measurements made on the enriched sample; and determining the likelihood that the man is the biological father of her fetus. The probability that the alleged father is the biological father of the fetus is determined to be 99.9999%, and the confidence of the paternity determination is calculated to be 99.99%. A report is printed out, or sent electronically to the pregnant woman's obstetrician, who transmits the determination to the woman. The woman, her husband, and the doctor sit down and discuss the report.

In an embodiment, the raw genetic material of the mother and the father is transformed by way of amplification to an amount of DNA that is similar in sequence, but larger in quantity. Then, by way of a genotyping method, the genotypic data that is encoded by nucleic acids is transformed into genetic measurements that may be stored physically and/or electronically on a memory device, such as those described above. The relevant algorithms that makeup the PARENTAL SUPPORT™ algorithm, relevant parts of which are discussed in detail herein, are translated into a computer program, using a programming language. Then, through the execution of the computer program on the computer hardware, instead of being physically encoded bits and bytes, arranged in a pattern that represents raw measurement data, they become transformed into a pattern that represents a high confidence determination of the paternity of the fetus. The details of this transformation will rely on the data itself and the computer language and hardware system used to execute the method described herein. Then, the data that is physically configured to represent a high quality paternity determination of the fetus is transformed into a report which may be sent to a health care practitioner. This transformation may be carried out using a printer or a computer display. The report may be a printed copy, on paper or other suitable medium, or else it may be electronic. In the case of an electronic report, it may be transmitted, it may be physically stored on a memory device at a location on the computer accessible by the health care practitioner; it also may be displayed on a screen so that it may be read. In the case of a screen display, the data may be transformed to a readable format by causing the physical transformation of pixels on the display device. The transformation may be accomplished by way of physically firing electrons at a phosphorescent screen, by way of altering an electric charge that physically changes the transparency of a specific set of pixels on a screen that may lie in front of a substrate that emits or absorbs photons. This transformation may be accomplished by way of changing the nanoscale orientation of the molecules in a liquid crystal, for example, from nematic to cholesteric or smectic phase, at a specific set of pixels. This transformation may be accomplished by way of an electric current causing photons to be emitted from a specific set of pixels made from a plurality of light emitting diodes arranged in a meaningful pattern. This transformation may be accomplished by any other way used to display information, such as a computer screen, or some other output device or way of transmitting information. The health care practitioner may then act on the report, such that the data in the report is transformed into an action. The action may be to continue or discontinue the pregnancy, in which case a gestating fetus is transformed into non-living fetus. Alternately, one may transform a set of genotypic measurements into a report that helps a physician treat his pregnant patient.

In some embodiments, the methods described herein can be used at a very early gestational age, for example as early as four week, as early as five weeks, as early as six weeks, as early as seven weeks, as early as eight weeks, as early as nine weeks, as early as ten weeks, as early as eleven weeks, and as early as twelve weeks.

Any of the embodiments disclosed herein may be implemented in digital electronic circuitry, integrated circuitry, specially designed ASICs (application-specific integrated circuits), computer hardware, firmware, software, or in combinations thereof. Apparatus of the presently disclosed embodiments can be implemented in a computer program product tangibly embodied in a machine-readable storage device for execution by a programmable processor; and method steps of the presently disclosed embodiments can be performed by a programmable processor executing a program of instructions to perform functions of the presently disclosed embodiments by operating on input data and generating output. The presently disclosed embodiments can be implemented advantageously in one or more computer programs that are executable and/or interpretable on a programmable system including at least one programmable processor, which may be special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device. Each computer program can be implemented in a high-level procedural or object-oriented programming language or in assembly or machine language if desired; and in any case, the language can be a compiled or interpreted language. A computer program may be deployed in any form, including as a stand-alone program, or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program may be deployed to be executed or interpreted on one computer or on multiple computers at one site, or distributed across multiple sites and interconnected by a communication network.

Computer readable storage media, as used herein, refers to physical or tangible storage (as opposed to signals) and includes without limitation volatile and non-volatile, removable and non-removable media implemented in any method or technology for the tangible storage of information such as computer-readable instructions, data structures, program modules or other data. Computer readable storage media includes, but is not limited to, RAM, ROM, EPROM, EEPROM, flash memory or other solid state memory technology, CD-ROM, DVD, or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other physical or material medium which can be used to tangibly store the desired information or data or instructions and which can be accessed by a computer or processor.

Targeted Enrichment and Sequencing

The use of a technique to enrich a sample of DNA at a set of target loci followed by sequencing as part of a method for non-invasive prenatal allele calling or ploidy calling may confer a number of unexpected advantages. In some embodiments of the present disclosure, the method involves measuring genetic data for use with an informatics based method, such as PARENTAL SUPPORT™ (PS). The ultimate outcome of some of the embodiments is the actionable genetic data of an embryo or a fetus. There are many methods that may be used to measure the genetic data of the individual and/or the related individuals as part of embodied methods. In an embodiment, a method for enriching the concentration of a set of targeted alleles is disclosed herein, the method comprising one or more of the following steps: targeted amplification of genetic material, addition of loci specific oligonucleotide probes, ligation of specified DNA strands, isolation of sets of desired DNA, removal of unwanted components of a reaction, detection of certain sequences of DNA by hybridization, and detection of the sequence of one or a plurality of strands of DNA by DNA sequencing methods. In some cases the DNA strands may refer to target genetic material, in some cases they may refer to primers, in some cases they may refer to synthesized sequences, or combinations thereof. These steps may be carried out in a number of different orders. Given the highly variable nature of molecular biology, it is generally not obvious which methods, and which combinations of steps, will perform poorly, well, or best in various situations.

For example, a universal amplification step of the DNA prior to targeted amplification may confer several advantages, such as removing the risk of bottlenecking and reducing allelic bias. The DNA may be mixed an oligonucleotide probe that can hybridize with two neighboring regions of the target sequence, one on either side. After hybridization, the ends of the probe may be connected by adding a polymerase, a means for ligation, and any necessary reagents to allow the circularization of the probe. After circularization, an exonuclease may be added to digest to non-circularized genetic material, followed by detection of the circularized probe. The DNA may be mixed with PCR primers that can hybridize with two neighboring regions of the target sequence, one on either side. After hybridization, the ends of the probe may be connected by adding a polymerase, a means for ligation, and any necessary reagents to complete PCR amplification. Amplified or unamplified DNA may be targeted by hybrid capture probes that target a set of loci; after hybridization, the probe may be localized and separated from the mixture to provide a mixture of DNA that is enriched in target sequences.

In some embodiments the detection of the target genetic material may be done in a multiplexed fashion. The number of genetic target sequences that may be run in parallel can range from one to ten, ten to one hundred, one hundred to one thousand, one thousand to ten thousand, ten thousand to one hundred thousand, one hundred thousand to one million, or one million to ten million. Note that the prior art includes disclosures of successful multiplexed PCR reactions involving pools of up to about 50 or 100 primers, and not more. Prior attempts to multiplex more than 100 primers per pool have resulted in significant problems with unwanted side reactions such as primer-dimer formation.

In some embodiments, this method may be used to genotype a single cell, a small number of cells, two to five cells, six to ten cells, ten to twenty cells, twenty to fifty cell, fifty to one hundred cells, one hundred to one thousand cells, or a small amount of extracellular DNA, for example from one to ten picograms, from ten to one hundred picograms, from one hundred picograms to one nanogram, from one to ten nanograms, from ten to one hundred nanograms, or from one hundred nanograms to one microgram.

The use of a method to target certain loci followed by sequencing as part of a method for allele calling or ploidy calling may confer a number of unexpected advantages. Some methods by which DNA may be targeted, or preferentially enriched, include using circularizing probes, linked inverted probes (LIPs, MIPs), capture by hybridization methods such as SURESELECT, and targeted PCR or ligation-mediated PCR amplification strategies.

In some embodiments, a method of the present disclosure involves measuring genetic data for use with an informatics based method, such as PARENTAL SUPPORT™ (PS). PARENTAL SUPPORT™ is an informatics based approach to manipulating genetic data, aspects of which are described herein. The ultimate outcome of some of the embodiments is the actionable genetic data of an embryo or a fetus followed by a clinical decision based on the actionable data. The algorithms behind the PS method take the measured genetic data of the target individual, often an embryo or fetus, and the measured genetic data from related individuals, and are able to increase the accuracy with which the genetic state of the target individual is known. In an embodiment, the measured genetic data is used in the context of making paternity determinations during prenatal genetic diagnosis. There are many methods that may be used to measure the genetic data of the individual and/or the related individuals in the aforementioned contexts. The different methods comprise a number of steps, those steps often involving amplification of genetic material, addition of oligonucleotide probes, ligation of specified DNA strands, isolation of sets of desired DNA, removal of unwanted components of a reaction, detection of certain sequences of DNA by hybridization, detection of the sequence of one or a plurality of strands of DNA by DNA sequencing methods. In some cases the DNA strands may refer to target genetic material, in some cases they may refer to primers, in some cases they may refer to synthesized sequences, or combinations thereof. These steps may be carried out in a number of different orders. Given the highly variable nature of molecular biology, it is generally not obvious which methods, and which combinations of steps, will perform poorly, well, or best in various situations.

Some embodiments of the present disclosure involve the use of “Linked Inverted Probes” (LIPs), which have been previously described in the literature. LIPs is a generic term meant to encompass technologies that involve the creation of a circular molecule of DNA, where the probes are designed to hybridize to targeted region of DNA on either side of a targeted allele, such that addition of appropriate polymerases and/or ligases, and the appropriate conditions, buffers and other reagents, will complete the complementary, inverted region of DNA across the targeted allele to create a circular loop of DNA that captures the information found in the targeted allele. LIPs may also be called pre-circularized probes, pre-circularizing probes, or circularizing probes. The LIPs probe may be a linear DNA molecule between 50 and 500 nucleotides in length, and in an embodiment between 70 and 100 nucleotides in length; in some embodiments, it may be longer or shorter than described herein. Others embodiments of the present disclosure involve different incarnations, of the LIPs technology, such as Padlock Probes and Molecular Inversion Probes (MIPs).

One method to target specific locations for sequencing is to synthesize probes in which the 3′ and 5′ ends of the probes anneal to target DNA at locations adjacent to and on either side of the targeted region, in an inverted manner, such that the addition of DNA polymerase and DNA ligase results in extension from the 3′ end, adding bases to single stranded probe that are complementary to the target molecule (gap-fill), followed by ligation of the new 3′ end to the 5′ end of the original probe resulting in a circular DNA molecule that can be subsequently isolated from background DNA. The probe ends are designed to flank the targeted region of interest. One aspect of this approach is commonly called MIPS and has been used in conjunction with array technologies to determine the nature of the sequence filled in.

Ligation-mediated PCR is method of PCR used to preferentially enrich a sample of DNA by amplifying one or a plurality of loci in a mixture of DNA, the method comprising: obtaining a set of primer pairs, where each primer in the pair contains a target specific sequence and a non-target sequence, where the target specific sequence is designed to anneal to a target region, one upstream and one downstream from the polymorphic site; polymerization of the DNA from the 3-prime end of upstream primer to the fill the single strand region between it and the 5-prime end of the downstream primer with nucleotides complementary to the target molecule; ligation of the last polymerized base of the upstream primer to the adjacent 5-prime base of the downstream primer; and amplification of only polymerized and ligated molecules using the non-target sequences contained at the 5-prime end of the upstream primer and the 3-prime end of the downstream primer. Pairs of primers to distinct targets may be mixed in the same reaction. The non-target sequences serve as universal sequences such that all pairs of primers that have been successfully polymerized and ligated may be amplified with a single pair of amplification primers.

In an embodiment, a sample of DNA may be preferentially enriched using a capture by hybridization approach. Some examples of commercial capture by hybridization technologies include AGILENT's SURESELECT and ILLUMINA's TRUSEQ. In capture by hybridization, a set of oligonucleotides that is complimentary or mostly complimentary to the desired targeted sequences is allowed to hybridize to a mixture of DNA, and then physically separated from the mixture. Once the desired sequences have hybridized to the targeting oligonucleotides, the effect of physically removing the targeting oligonucleotides is to also remove the targeted sequences. Once the hybridized oligos are removed, they can be heated to above their melting temperature and they can be amplified. Some ways to physically remove the targeting oligonucleotides is by covalently bonding the targeting oligos to a solid support, for example a magnetic bead, or a chip. Another way to physically remove the targeting oligonucleotides is by covalently bonding them to a molecular moiety with a strong affinity for another molecular moiety. An example of such a molecular pair is biotin and streptavidin, such as is used in SURESELECT. Thus that targeted sequences could be covalently attached to a biotin molecule, and after hybridization, a solid support with streptavidin affixed can be used to pull down the biotinylated oligonucleotides, to which are hybridized to the targeted sequences.

In some embodiments, PCR can be used to target specific locations of the genome. In plasma samples, the original DNA is highly fragmented (typically less than 500 bp, with an average length less than 200 bp). In PCR, both forward and reverse primers must anneal to the same fragment to enable amplification. Therefore, if the fragments are short, the PCR assays must amplify relatively short regions as well. PCR assay can be generated in large numbers, however, the interactions between different PCR assays makes it difficult to multiplex them beyond about one hundred assays. Various complex molecular approaches can be used to increase the level of multiplexing, but it may still be limited to fewer than 100, perhaps 200, or possibly 500 assays per reaction. Samples with large quantities of DNA can be split among multiple sub-reactions and then recombined before sequencing. For samples where either the overall sample or some subpopulation of DNA molecules is limited, splitting the sample would introduce statistical noise. In an embodiment, a small or limited quantity of DNA may refer to an amount below 10 pg, between 10 and 100 pg, between 100 pg and 1 ng, between 1 and 10 ng, or between 10 and 100 ng. Note that while this method is particularly useful on small amounts of DNA where other methods that involve splitting into multiple pools can cause significant problems related to introduced stochastic noise, this method still provides the benefit of minimizing bias when it is run on samples of any quantity of DNA. In these situations a universal pre-amplification step may be used to increase the overall sample quantity. Ideally, this pre-amplification step should not appreciably alter the allelic distributions.

In general, to perform targeted sequencing of multiple (n) targets of a sample (greater than 50, greater than 100, greater than 500, or greater than 1,000), one can split the sample into a number of parallel reactions that amplify one or a smaller number of individual targets. This has been performed in PCR multiwell plates or can be done in commercial platforms such as the FLUIDIGM ACCESS ARRAY (48 reactions per sample in microfluidic chips) or DROPLET PCR by RAIN DANCE TECHNOLOGY (100s to a few thousands of targets). Unfortunately, these split-and-pool methods are problematic for samples with a limited amount of DNA, as there is often not enough copies of the genome to ensure that there is one copy of each region of the genome in each well. This is an especially severe problem when polymorphic loci are targeted, and the relative proportions of the alleles at the polymorphic loci are needed, as the stochastic noise introduced by the splitting and pooling will cause very poorly accurate measurements of the proportions of the alleles that were present in the original sample of DNA. Described here is a method to effectively and efficiently amplify many PCR reactions that is applicable to cases where only a limited amount of DNA is available. In an embodiment, the method may be applied for analysis of single cells, body fluids, mixtures of DNA such as the free floating DNA found in maternal plasma, biopsies, environmental and/or forensic samples.

In an embodiment, the targeted sequencing may involve one, a plurality, or all of the following steps. a) Generate and universally amplify a library with adaptor sequences on both ends of DNA fragments. b) Divide into multiple reactions after library amplification. c) Perform about 100-plex, about 1000-plex, or about 10,000-plex amplification of selected targets using one target specific “Forward” primer per target and one tag specific primer. d) Perform a second amplification from this product using “Reverse” target specific primers and one (or more) primer specific to a universal tag that was introduced as part of the target specific forward primers in the first round. e) Divide the product into multiple aliquots and amplify subpools of targets in individual reactions (for example, 50 to 500-plex, though this can be used all the way down to singleplex. f) Pool products of parallel subpools reactions. During these amplifications primers may carry sequencing compatible tags (partial or full length) such that the products can be sequenced.

In an embodiment, it is possible to mitigate potential losses in subsequent steps by amplifying all or a fraction of the original cell free DNA (cfDNA) sample. Various methods are available to amplify all of the genetic material in a sample, increasing the amount available for downstream procedures. In an embodiment, ligation mediated PCR (LM-PCR) DNA fragments are amplified by PCR after ligation of either one distinct adaptors, two distinct adapters, or many distinct adaptors. In an embodiment, multiple displacement amplification (MDA) phi-29 polymerase is used to amplify all DNA isothermally. In DOP-PCR and variations, random priming is used to amplify the original material DNA. Each method has certain characteristics such as uniformity of amplification across all represented regions of the genome, efficiency of capture and amplification of original DNA, and amplification performance as a function of the length of the fragment.

Traditional PCR assay design results in significant losses of distinct fetal molecules, but losses can be greatly reduced by designing very short PCR assays, termed mini-PCR assays. Fetal cfDNA in maternal serum is highly fragmented and the fragment sizes are distributed in approximately a Gaussian fashion with a mean of about 160 bp, a standard deviation of about 15 bp, a minimum size of about 100 bp, and a maximum size of about 220 bp. The distribution of fragment start and end positions with respect to the targeted polymorphisms, while not necessarily random, vary widely among individual targets and among all targets collectively and the polymorphic site of one particular target locus may occupy any position from the start to the end among the various fragments originating from that locus. Note that the term mini-PCR may equally well refer to normal PCR with no additional restrictions or limitations.

During PCR, amplification will only occur from template DNA fragments comprising both forward and reverse primer sites. Because fetal cfDNA fragments are short, the likelihood of both primer sites being present the likelihood of a fetal fragment of length L comprising both the forward and reverse primers sites is ratio of the length of the amplicon to the length of the fragment. Under ideal conditions, assays in which the amplicon is 45, 50, 55, 60, 65, or 70 bp will successfully amplify from about 72%, 69%, 66%, 63%, 59%, or 56%, respectively, of available template fragment molecules. The amplicon length is the distance between the 5-prime ends of the forward and reverse priming sites. Amplicon length that is shorter than typically used by those known in the art may result in more efficient measurements of the desired polymorphic loci by only requiring short sequence reads. In an embodiment, a substantial fraction of the amplicons should be less than 100 bp, less than 90 bp, less than 80 bp, less than 70 bp, less than 65 bp, less than 60 bp, less than 55 bp, less than 50 bp, or less than 45 bp.

Note that in methods known in the prior art, short assays such as those described herein are usually avoided because they are not required and they impose considerable constraint on primer design by limiting primer length, annealing characteristics, and the distance between the forward and reverse primer.

Multiplex PCR may involve a single round of PCR in which all targets are amplified or it may involve one round of PCR followed by one or more rounds of nested PCR or some variant of nested PCR. Nested PCR consists of a subsequent round or rounds of PCR amplification using one or more new primers that bind internally, by at least one base pair, to the primers used in a previous round. Nested PCR reduces the number of spurious amplification targets by amplifying, in subsequent reactions, only those amplification products from the previous one that have the correct internal sequence. Reducing spurious amplification targets improves the number of useful measurements that can be obtained, especially in sequencing. Nested PCR typically entails designing primers completely internal to the previous primer binding sites, necessarily increasing the minimum DNA segment size required for amplification. For samples such as maternal plasma cfDNA, in which the DNA is highly fragmented, the larger assay size reduces the number of distinct cfDNA molecules from which a measurement can be obtained. In an embodiment, to offset this effect, one may use a partial nesting approach where one or both of the second round primers overlap the first binding sites extending internally some number of bases to achieve additional specificity while minimally increasing in the total assay size.

In an embodiment, a multiplex pool of PCR assays are designed to amplify potentially heterozygous SNP or other polymorphic or non-polymorphic loci on one or more chromosomes and these assays are used in a single reaction to amplify DNA. The number of PCR assays may be between 50 and 200 PCR assays, between 200 and 1,000 PCR assays, between 1,000 and 5,000 PCR assays, or between 5,000 and 20,000 PCR assays (50 to 200-plex, 200 to 1,000-plex, 1,000 to 5,000-plex, 5,000 to 20,000-plex, more than 20,000-plex respectively).

In an embodiment, a 100-plex to -500 plex, 500-plex to 1,000-plex, 1,000-plex to 2,000-plex, 2,000-plex to 5,000-plex, 5,000-plex to 10,000-plex, 10,000-plex to 20,000-plex, 20,000-plex to 50,000-plex, or 50,000-plex to 100,000-plex PCR assay pool is created such that forward and reverse primers have tails corresponding to the required forward and reverse sequences required by a high throughput sequencing instrument such as the HISEQ, GAIIX, or MISEQ available from ILLUMINA. In addition, included 5-prime to the sequencing tails is an additional sequence that can be used as a priming site in a subsequent PCR to add nucleotide barcode sequences to the amplicons, enabling multiplex sequencing of multiple samples in a single lane of the high throughput sequencing instrument. In an embodiment, a 10,000-plex PCR assay pool is created such that reverse primers have tails corresponding to the required reverse sequences required by a high throughput sequencing instrument. After amplification with the first 10,000-plex assay, a subsequent PCR amplification may be performed using a another 10,000-plex pool having partly nested forward primers (e.g. 6-bases nested) for all targets and a reverse primer corresponding to the reverse sequencing tail included in the first round. This subsequent round of partly nested amplification with just one target specific primer and a universal primer limits the required size of the assay, reducing sampling noise, but greatly reduces the number of spurious amplicons. The sequencing tags can be added to appended ligation adaptors and/or as part of PCR probes, such that the tag is part of the final amplicon.

Fetal fraction affects performance of the test; it is more difficult to determine the correct paternity on samples with a lower fetal fraction. There are a number of ways to enrich the fetal fraction of the DNA found in maternal plasma. Fetal fraction can be increased by the previously described LM-PCR method already discussed as well as by a targeted removal of long maternal fragments. In embodiment, the longer fragments are removed using size selection techniques. In an embodiment, prior to multiplex PCR amplification of the target loci, an additional multiplex PCR reaction may be carried out to selectively remove long and largely maternal fragments corresponding to the loci targeted in the subsequent multiplex PCR. Additional primers are designed to anneal a site a greater distance from the polymorphism than is expected to be present among cell free fetal DNA fragments. These primers may be used in a one cycle multiplex PCR reaction prior to multiplex PCR of the target polymorphic loci. These distal primers are tagged with a molecule or moiety that can allow selective recognition of the tagged pieces of DNA. In an embodiment, these molecules of DNA may be covalently modified with a biotin molecule that allows removal of newly formed double stranded DNA comprising these primers after one cycle of PCR. Double stranded DNA formed during that first round is likely maternal in origin. Removal of the hybrid material may be accomplish by the used of magnetic streptavidin beads. There are other methods of tagging that may work equally well. In an embodiment, size selection methods may be used to enrich the sample for shorter strands of DNA; for example those less than about 800 bp, less than about 500 bp, or less than about 300 bp. Amplification of short fragments can then proceed as usual.

In some embodiments, the target DNA may originate from single cells, from samples of DNA consisting of less than one copy of the target genome, from low amounts of DNA, from DNA from mixed origin (e.g. pregnancy plasma: placental and maternal DNA; cancer patient plasma and tumors: mix between healthy and cancer DNA, transplantation etc), from other body fluids, from cell cultures, from culture supernatants, from forensic samples of DNA, from ancient samples of DNA (e.g. insects trapped in amber), from other samples of DNA, and combinations thereof.

In some embodiments, the methods described herein may be used to amplify and/or detect SNPs, single tandem repeats (STRs), copy number, nucleotide methylation, mRNA levels, other types of RNA expression levels, other genetic and/or epigenetic features. The mini-PCR methods described herein may be used along with next-generation sequencing; it may be used with other downstream methods such as microarrays, counting by digital PCR, real-time PCR, Mass-spectrometry analysis etc.

In some embodiment, the mini-PCR amplification methods described herein may be used as part of a method for accurate quantification of minority populations. It may be used for absolute quantification using spike calibrators. It may be used for mutation/minor allele quantification through very deep sequencing, and may be run in a highly multiplexed fashion. It may be used for standard paternity and identity testing of relatives or ancestors, in human, animals, plants or other creatures. It may be used for forensic testing. It may be used for rapid genotyping and copy number analysis (CN), on any kind of material, e.g. amniotic fluid and CVS, sperm, product of conception (POC). It may be used for single cell analysis, such as genotyping on samples biopsied from embryos. It may be used for rapid embryo analysis (within less than one, one, or two days of biopsy) by targeted sequencing using min-PCR.

Highly multiplexed PCR can often result in the production of a very high proportion of product DNA resulting from unproductive side reactions such as primer dimer formation. In an embodiment, the particular primers that are most likely to cause unproductive side reactions may be removed from the primer library to give a primer library that will result in a greater proportion of amplified DNA that maps to the genome. The step of removing problematic primers, that is, those primers that are particularly likely to firm dimers has unexpectedly enabled extremely high PCR multiplexing levels for subsequent analysis by sequencing. In systems such as sequencing, where performance significantly degrades by primer dimers and/or other mischief products, greater than 10, greater than 50, and greater than 100 times higher multiplexing than other described multiplexing has been achieved. Note this is opposed to probe based detection methods, e.g. microarrays, TaqMan, PCR etc. where an excess of primer dimers will not affect the outcome appreciably. Also note that the general belief in the art is that multiplexing PCR for sequencing is limited to about 100 assays in the same well.

There are a number of ways to choose primers for a library where the amount of non-mapping primer-dimer or other primer mischief products are minimized. Empirical data indicate that a small number of ‘bad’ primers are responsible for a large amount of non-mapping primer dimer side reactions. Removing these ‘bad’ primers can increase the percent of sequence reads that map to targeted loci. One way to identify the ‘bad’ primers is to look at the sequencing data of DNA that was amplified by targeted amplification; those primer dimers that are seen with greatest frequency can be removed to give a primer library that is significantly less likely to result in side product DNA that does not map to the genome. There are also publicly available programs that can calculate the binding energy of various primer combinations, and removing those with the highest binding energy will also give a primer library that is significantly less likely to result in side product DNA that does not map to the genome.

Note that there are other methods for determining which PCR probes are likely to form dimers. In an embodiment, analysis of a pool of DNA that has been amplified using a non-optimized set of primers may be sufficient to determine problematic primers. For example, analysis may be done using sequencing, and those dimers which are present in the greatest number are determined to be those most likely to form dimers, and may be removed.

To select target locations, one may start with a pool of alleged primer pair designs and create a thermodynamic model of potentially adverse interactions between primer pairs, and then use the model to eliminate designs that are incompatible with other the designs in the pool.

There are many workflows that are possible when conducting PCR; some workflows typical to the methods disclosed herein are described. The steps outlined herein are not meant to exclude other possible steps nor does it imply that any of the steps described herein are required for the method to work properly. A large number of parameter variations or other modifications are known in the literature, and may be made without affecting the essence of the invention. One particular generalized workflow is given below followed by a number of possible variants. The variants typically refer to possible secondary PCR reactions, for example different types of nesting that may be done (step 3). It is important to note that variants may be done at different times, or in different orders than explicitly described herein.

1. The DNA in the sample may have ligation adapters, often referred to as library tags or ligation adaptor tags (LTs), appended, where the ligation adapters contain a universal priming sequence, followed by a universal amplification. In an embodiment, this may be done using a standard protocol designed to create sequencing libraries after fragmentation. In an embodiment, the DNA sample can be blunt ended, and then an A can be added at the 3′ end. A Y-adaptor with a T-overhang can be added and ligated. In some embodiments, other sticky ends can be used other than an A or T overhang. In some embodiments, other adaptors can be added, for example looped ligation adaptors. In some embodiments, the adaptors may have tag designed for PCR amplification. 2. Specific Target Amplification (STA): Pre-amplification of hundreds to thousands to tens of thousands and even hundreds of thousands of targets may be multiplexed in one reaction. STA is typically run from 10 to 30 cycles, though it may be run from 5 to 40 cycles, from 2 to 50 cycles, and even from 1 to 100 cycles. Primers may be tailed, for example for a simpler workflow or to avoid sequencing of a large proportion of dimers. Note that typically, dimers of both primers carrying the same tag will not be amplified or sequenced efficiently. In some embodiments, between 1 and 10 cycles of PCR may be carried out; in some embodiments between 10 and 20 cycles of PCR may be carried out; in some embodiments between 20 and 30 cycles of PCR may be carried out; in some embodiments between 30 and 40 cycles of PCR may be carried out; in some embodiments more than 40 cycles of PCR may be carried out. The amplification may be a linear amplification. The number of PCR cycles may be optimized to result in an optimal depth of read (DOR) profile. Different DOR profiles may be desirable for different purposes. In some embodiments, a more even distribution of reads between all assays is desirable; if the DOR is too small for some assays, the stochastic noise can be too high for the data to be too useful, while if the depth of read is too high, the marginal usefulness of each additional read is relatively small.

Primer tails may improve the detection of fragmented DNA from universally tagged libraries. If the library tag and the primer-tails contain a homologous sequence, hybridization can be improved (for example, melting temperature (T_(M)) is lowered) and primers can be extended if only a portion of the primer target sequence is in the sample DNA fragment. In some embodiments, 13 or more target specific base pairs may be used. In some embodiments, 10 to 12 target specific base pairs may be used. In some embodiments, 8 to 9 target specific base pairs may be used. In some embodiments, 6 to 7 target specific base pairs may be used. In some embodiments, STA may be performed on pre-amplified DNA, e.g. MDA, RCA, other whole genome amplifications, or adaptor-mediated universal PCR. In some embodiments, STA may be performed on samples that are enriched or depleted of certain sequences and populations, e.g. by size selection, target capture, directed degradation.

3. In some embodiments, it is possible to perform secondary multiplex PCRs or primer extension reactions to increase specificity and reduce undesirable products. For example, full nesting, semi-nesting, hemi-nesting, and/or subdividing into parallel reactions of smaller assay pools are all techniques that may be used to increase specificity. Experiments have shown that splitting a sample into three 400-plex reactions resulted in product DNA with greater specificity than one 1,200-plex reaction with exactly the same primers. Similarly, experiments have shown that splitting a sample into four 2,400-plex reactions resulted in product DNA with greater specificity than one 9,600-plex reaction with exactly the same primers. In an embodiment, it is possible to use target-specific and tag specific primers of the same and opposing directionality. 4. In some embodiments, it is possible to amplify a DNA sample (dilution, purified or otherwise) produced by an STA reaction using tag-specific primers and “universal amplification”, i.e. to amplify many or all pre-amplified and tagged targets. Primers may contain additional functional sequences, e.g. barcodes, or a full adaptor sequence necessary for sequencing on a high throughput sequencing platform.

These methods may be used for analysis of any sample of DNA, and are especially useful when the sample of DNA is particularly small, or when it is a sample of DNA where the DNA originates from more than one individual, such as in the case of maternal plasma. These methods may be used on DNA samples such as a single or small number of cells, genomic DNA, plasma DNA, amplified plasma libraries, amplified apoptotic supernatant libraries, or other samples of mixed DNA. In an embodiment, these methods may be used in the case where cells of different genetic constitution may be present in a single individual, such as with cancer or transplants.

In some embodiments, the multiplex PCR amplification may involve using various types of nesting protocol, for example: semi-nested mini-PCR, fully nested mini-PCR, heminested mini-PCR, triply hemi-nested mini-PCR, one-sided nested mini-PCR, one-sided mini-PCR or reverse semi-nested mini-PCR.

Diagnostic Box

In an embodiment, the present disclosure comprises a diagnostic box that is capable of partly or completely carrying out aspects of the methods described in this disclosure. In an embodiment, the diagnostic box may be located at a physician's office, a hospital laboratory, or any suitable location reasonably proximal to the point of patient care. The box may be able to run the aspects of the method in a wholly automated fashion, or the box may require one or a number of steps to be completed manually by a technician. In an embodiment, the box may be able to analyze the genotypic data measured on the maternal plasma. In an embodiment, the box may be linked to means to transmit the genotypic data measured using the diagnostic box to an external computation facility which may then analyze the genotypic data, and possibly also generate a report. The diagnostic box may include a robotic unit that is capable of transferring aqueous or liquid samples from one container to another. It may comprise a number of reagents, both solid and liquid. It may comprise a high throughput sequencer. It may comprise a computer.

Primer Kit

In some embodiments, a kit may be formulated that comprises a plurality of primers designed to achieve the methods described in this disclosure. The primers may be outer forward and reverse primers, inner forward and reverse primers as disclosed herein; they could be primers that have been designed to have low binding affinity to other primers in the kit as disclosed in the section on primer design; they could be hybrid capture probes or pre-circularized probes as described in the relevant sections, or some combination thereof. In an embodiment, a kit may be formulated for determining a ploidy status of a target chromosome in a gestating fetus designed to be used with the methods disclosed herein, the kit comprising a plurality of inner forward primers and optionally the plurality of inner reverse primers, and optionally outer forward primers and outer reverse primers, where each of the primers is designed to hybridize to the region of DNA immediately upstream and/or downstream from one of the polymorphic sites on the target chromosome, and optionally additional chromosomes. In an embodiment, the primer kit may be used in combination with the diagnostic box described elsewhere in this document.

Maximum Likelihood Estimates

Many methods known in the art for detecting the presence or absence of a phenotype or genotype, for example, a chromosomal abnormality, a medical condition or a paternity relationship involve the use of a single hypothesis rejection test, where a metric that is directly related to or correlated with the condition is measured, and if the metric is on one side of a given threshold, the condition is determined to be present, while if the metric falls on the other side of the threshold, the condition is determined to be absent. A single-hypothesis rejection test only looks at the null distribution when deciding between the null and alternate hypotheses. Without taking into account the alternate distribution, one cannot estimate the likelihood of each hypothesis given the observed data and therefore cannot calculate a confidence on the call. Hence with a single-hypothesis rejection test, one gets a yes or no answer without an estimate of the confidence associated with the specific case.

In some embodiments, the method disclosed herein is able to detect the presence or absence of phenotype or genotype, for example, a chromosomal abnormality, a medical condition or a paternity relationship, using a maximum likelihood method. This is a substantial improvement over a method using a single hypothesis rejection technique as the threshold for calling absence or presence of the condition can be adjusted as appropriate for each case. This is particularly relevant for diagnostic techniques that aim to determine the paternity of a gestating fetus from genetic data available from the mixture of fetal and maternal DNA present in the free floating DNA found in maternal plasma. The maximum likelihood estimation method may use the allelic distributions associated with each hypothesis to estimate the likelihood of the data conditioned on each hypothesis. These conditional probabilities can then be converted to a hypothesis call and confidence. Similarly, maximum a posteriori estimation method uses the same conditional probabilities as the maximum likelihood estimate, but also incorporates population priors when choosing the best hypothesis and determining confidence.

Therefore, the use of a maximum likelihood estimate (MLE) technique, or the closely related maximum a posteriori (MAP) technique give two advantages, first it increases the chance of a correct call, and it also allows a confidence to be calculated for each call. In an embodiment, selecting the paternity call corresponding to the hypothesis with the greatest probability is carried out using maximum likelihood estimates or maximum a posteriori estimates. In an embodiment, a method is disclosed for determining the paternity of a gestating fetus that involves taking any method currently known in the art that uses a single hypothesis rejection technique and reformulating it such that it uses a MLE or MAP technique

A Method for Paternity Determination

The ffDNA is typically present at low fraction in a mixture with maternal DNA. In some embodiments, the mother has known genotype or the maternal genotype can be measured or inferred. Typically, the fraction of fetal DNA found in maternal plasma is between 2 and 20%, although in different conditions this percentage can range from about 0.01% to about 50%. In an embodiment, a microarray or other technique that gives intensity data on a per allele basis can be used to measure the maternal plasma. In an embodiment, sequencing can be used to measure the DNA contained in the maternal plasma. In these cases the allele intensity measurement or sequence read count at a particular allele is a sum of the maternal and fetal signals. Assuming that the mixture ratio of child to mother DNA is r to 1, the relative number of alleles at a locus consists of 2 alleles from the mother and 2r alleles from the child. In some embodiment, the loci comprise single nucleotide polymorphisms. Table 1 shows the relative number of each allele in the mixture for a selection of informative parent contexts.

TABLE 1 Number of alleles by context parent context A in mixture B in mixture AA|AA 2 + 2r 0 AA|BB 2 + r  r AB|AA 1 + 2r 1 AA|AB 2 + 2r or 2 + r 0 or r BB|AA r 2 + r  BB|BB 0 2 + 2r

Note that the choice of the above four contexts as being informative is not meant to be inclusive of all contexts that may be informative. Any combination of contexts may be used, and there is a significant amount of information that may be found in the genotypic measurements of any context.

Even in the presence of significant allele dropout rate from the child, there may be a clear distinction between signal where an allele is present and signal where an allele is not present. For example, consider the A allele measurements from SNPs in context pairs BB|BB and BB|AA and the B allele measurements from SNPs in from context pairs AA|AA and AA|BB. In each case, there should be no signal present in the first context and there should be signal present in the second context, wherever the child's alleles have not dropped out. However, if the alleged father is incorrect, there will sometimes be signal present in both contexts. Thus, the distribution of SNP measurements should be different, depending on whether the alleged father is correct or not.

The difference will typically be more observable at the high-signal end of the distribution, because these will be the SNPs where there is higher likelihood of having DNA contributions from the child. This difference can be observed by comparing high percentile values of the distributions of SNP measurements. Examples of possible percentile methods are nearest rank, linear interpolation between closest ranks, and weighted percentile.

For example, define X₁ as the set of A allele SNP measurements in context BB|BB and X₂ as the set of A allele measurements in context BB|AA, from all chromosomes. If the alleged father is correct, then the 99^(th) percentile value of X₁ will be significantly less than the 99^(th) percentile value of X₂. If the alleged father is incorrect, the 99^(th) percentile values of the two distributions will be closer together. Note than any percentile may be used equally well, for example, 95^(th) percentile, 90^(th) percentile, 85^(th) percentile or 80^(th) percentile. In an embodiment, for a particular measurement channel, X1 can be defined as the measurements from the context with no signal and X2 can be defined as the measurements from the context where the mother and father are both homozygous, and only the father alleles provide a signal (inherited through the child).

Define p₁ as the 99th (or 95^(th), 90^(th) etc.) percentile of the X₁ data and p₂ as the 99^(th) (or 95^(th), 90^(th) etc.) percentile of the X₂ data. Define the test statistic t as p₁/p₂. Note that other functions of p1 and p2 that demonstrate the difference in values may be used equally well. The value of t will vary depending on the amount of child DNA present in the sample, which is not known. Therefore, classification thresholds for t can not be calculated a priori.

The test statistic t for a single sample can be compared to a distribution generated from the genotypes of many individuals who are known not to be the father, using the following procedure. Assume that the genotypes from a large set of unrelated individuals are available.

-   -   1. For each unrelated individual, assume that it is the father         and calculate the value of the test statistic t.     -   2. Let T_(u) be the set of t measurements from the unrelated         males. Fit a distribution to T_(u). This is the distribution of         t for the particular sample, under the null hypothesis. The null         hypothesis is “genotypes do not come from father of child         present in sample”. The distribution Pu(t) could be a maximum         likelihood fit or method of moments fit to a known distribution         e.g. a Gaussian distribution, a kernel density fit using a         kernel function e.g. Gaussian, box, triangle etc, or any other         appropriate distribution fitting method.     -   3. Consider the genotypes of the alleged father and calculate         the corresponding test statistic t_(c).     -   4. The true father is expected to result in a smaller value of t         than an unrelated individual. The probability of an unrelated         father producing t_(c) or a more outlying value is the         cumulative density function of P_(u) evaluated at t_(c). Thus,         the p-value p for rejecting the null hypothesis is given by the         following:

p=∫ ₀ ^(t) ^(c) P _(u)(t)dt

If p falls below a significance threshold α then the hypothesis that the father alleged is an unrelated individual can be rejected with significance α. If p is greater than α then the null hypothesis cannot be rejected, and the alleged father may be unrelated to the child present in the sample. The significance threshold α defines the probability that an unrelated individual could be classified as the correct father. For example, with a threshold of a equal to 0.01, one percent of unrelated individuals could be identified as the correct father.

Various methods can be considered for combination of data from the A and B allele channels. For example, two simple methods are to require that the p-value from all channels be below a threshold, or to require that the p-value from any channel be below the threshold.

In some embodiments, the paternity testing method assumes that child DNA is present at sufficient concentration to distinguish between SNPs that have or do not have signal from the child. In the absence of sufficient child DNA concentration, this method may report “incorrect father” because expected paternally inherited alleles are not measured in the maternal plasma. In an embodiment, a method is described that can confirm the presence of sufficient child DNA before applying the paternity test. The child presence confirmation is based on calculation of a test statistic that is proportional to the child DNA concentration, but does not require father genotypes. If the test statistic is above the required threshold, then the concentration of child DNA is sufficient to perform the paternity test.

Consider the set of SNPs (from all chromosomes combined) where the mother genotype is AA and the B channel is measured. A signal is expected only on the subset of SNPs where the child genotype contains a B, but these SNPs cannot be identified a priori without the father genotypes, which are not available. Instead, consider the SNP populations frequencies {f_(i)} where f_(i) is the sample mean number of Bs in the genotype of SNP i, based on a large sample population. Note that most SNPs where the mother genotype is AA will have f_(i) less than 0.5, but the distribution of f_(i) on these SNPs extends almost to one. Consider two sets of SNPs, S₁ and S₂, where S₁={i: f_(i)<T_(L)} and S2={i: f_(i)>T_(H)}. The thresholds T_(L) and T_(H) are set so that very few SNPs in S1 are expected to have a B and many SNPs in S2 are expected to have a B, and each set has sufficient population. In one embodiment, the algorithm uses T_(L)=0.05 and T_(H)=0.7, while other values for TL and TH might work equally well or better. Let y_(i) be the B channel measurement from SNP i, Y₁={y_(i): iϵS₁}, and Y₂={y_(i): iϵS₂}. The distributions of Y₁ and Y₂ will be very similar because most SNPs in both distributions will have no signal. However, some non-trivial number of SNPs in S₂ are expected to have child signal and very few SNPs in S₁ are expected to have child signal. Therefore, the tail of Y₂ should extend to higher intensity than the tail of Y₁. Let p_(t) be a percentile close to 1, for example, the 99th percentile. In the presence of sufficient child DNA, the p_(t) percentile of Y₂ should be significantly higher than the p_(t) percentile of Y₁. Thus, the test statistic s can be defined as follows.

s=percentile(Y ₂ ;p _(t))−percentile(Y ₁ ;p _(t))

In one embodiment, the test statistic may be normalized by a variety of methods to attempt to account for amplification differences between arrays. In one embodiment, the normalization could be done on a per chromosome basis. In one embodiment, the normalization could be done on a per array basis. In one embodiment, the normalization could be done on a per sequencing run basis.

The following calculation shows how the thresholds T_(L) and T_(H) are able to distinguish the effect of child DNA in particular maternal sample, based on approximate numbers of SNPs and dropout rates. Table 3 shows some data from

In one embodiment, the method involves the following assumptions:

-   -   Population frequencies are calculated from a large population         data set, for example, more than 500 individuals, more than         1,000 individuals, more than 5,000 individuals or more than         20,000 individuals, and the number of SNPs in each context comes         from an example mother and father.     -   There are no SNPs where mother is AA and father is BB (in         reality, these are approximately 8 percent of mother AA SNPs)     -   Half of SNPs where father has B result in child B.     -   Child dropout rate is 90 percent.     -   Measurements with child signal will be higher than measurements         without child signal

Table 3 shows some data from a particular paternity case using the disclosed method with the above parameters. The 98th percentile measurement from S₁ is not expected to include any SNPs with child signal present. The 98th percentile measurement from S2 is expected to include about 50 SNPs with child signal present. The difference between the two should reflect the amount of child signal.

TABLE 3 Data pertaining to paternity determination num SNPs average f_(i) num SNPs num SNPs fraction of set, set definition in set in set father B child signal child signal S₁ f_(i) < 0.05 13300 0.012 171 9 0.0007 S₂ f_(i) > 0.7 3000 0.79 2370 119 0.039

FIG. 1 shows the distribution of allele intensity data for contexts AA|AA and AA|BB from a maternal plasma sample collected at 38 weeks. The B allele is measured. Note that the AA|BB distribution extends significantly higher than the AA|AA distribution, showing that the B allele (which is only present in the child's genome) is present in the AA|BB context but not the AA|AA context. FIG. 2 comes from the same maternal plasma sample as FIG. 1, and shows the distribution of the test statistic t for allele B using the genotypes from 200 unrelated individuals. Two distributions are shown (the two curves): the maximum likelihood Gaussian fit and a kernel distribution. The value t_(c) for the biological father is marked with a star. The p-value is less than 10⁻⁷ for the null hypothesis that the alleged father is unrelated to the child.

Table 4 presents results from 8 maternal blood samples at varying stages of pregnancy. A p-value is calculated based on the data measured from each channel (A allele and B allele) for the correct father. If both channel p-values are required to be below 0.01, then two samples are classified incorrectly. If only one channel is required to pass the threshold, then all samples are correctly classified. Any number of metrics and thresholds may be used for confirming or excluding parentage. For example, one could use a cut off p-value of 0.02, 0.005, 0.001, or 0.00001; similarly, one could demand that one or both channel p-values are below a given threshold, or one could have two different thresholds for the different channel p-values.

TABLE 4 P-values for two channels for eight paternity determinations. Weeks pregnancy p-value (Y) p-value (X) 11 2.3 × 10⁻⁷ <10⁻⁷ 16 0.013 <10⁻⁷ 17 <10⁻⁷ <10⁻⁷ 17 <10⁻⁷ 0.0002 20 <10⁻⁷ <10⁻⁷ 28 0.14 0.0048 38 <10⁻⁷ <10⁻⁷ 38 <10⁻⁷ <10⁻⁷

Table 4 shows the P-values for the null hypothesis that the correct father is an unrelated individual. Each row corresponds to a different maternal blood sample, and the corresponding paternal genetic sample. Genetic measurements made on 200 unrelated males were used as a control. The curve in FIG. 3 shows the distribution of intensity ratios for 200 unrelated males, and the star represents the intensity ratio for the biological father. This data is taken from a case where the blood was drawn from a mother who was 11 weeks pregnant.

FIG. 4 shows the cumulative distribution frequency (cdf) curves for the correlation ratio between the fetal genotypic measurements and the parental genotypic measurements for three cases: (1) where both the pregnant mother and the alleged father are the biological parents of the fetus (“correct”, rightmost curve), (2) where the pregnant mother is the biological mother of the fetus, but the alleged father is not the biological father of the fetus (“one wrong”, middle curve) and (3) where neither the pregnant mother nor the alleged father are the biological parents of the fetus (“two wrong”, leftmost curve). The cdf curves are the correlation ratio between genotypic data of the embryo, calculated from data measured on a single cell, and the genotypic data of the assumed parents when zero, one or two of the assumed parents are actually the genetic parents of the fetus. Note that the labels for “correct” and “two wrong” are reversed. FIG. 5 shows histograms for the same three cases. Note that this histogram is made up of more than 1000 cases where one or both parents are incorrect. The histogram of correlation rate measured between the genotypic data of the fetus, as measured on a single cell, and the genotypic data of the assumed parents when zero, one or two of the assumed parents are actually the genetic parents of the fetus.

Thirty five paternity results are shown in FIG. 6 using the instant method for paternity testing. They were run on samples collected from pregnant women with gestational ages ranging from 9 to 40 weeks. The red curve on the right represents a normalized Gaussian distribution of the paternity testing statistic for 800 unrelated males. The distribution of unrelated males is different for each case; a normalized distribution is used here for visualization purposes.

The blue bars represent the normalized test statistic for the correct (suspected) father. It is clear that the correct fathers are clearly separated from the unrelated males. Note that the normalized test statistic crudely approximates standard deviations, therefore, “−5” on the graph below is about 5 standard deviations from the mean. Thus all assumed correct fathers in this cohort have been confirmed as the correct fathers with a significance of at least 99.9999%.

In one embodiment of the invention, the knowledge of the parental haplotypes could be used to increase the accuracy of the test. For example, if the two haplotypes of the father are known for a given segment of a chromosome, then the knowledge of which SNPs are present for cases where there is no drop out can be used to determine which SNPs should be expected for those cases where there may be dropout. For example, imagine a set of three SNPs that are linked, that is, they are located close together on the same chromosome, and where the contexts of the mother and the alleged father are: AA|AB, AA|AB, AA|BA. Note that when the genotype of a parent is phased, then AB≠BA, since the first of the two letters represents the alleles on the first haplotype, and the second letter represents the alleles on the second haplotype. Now imagine that for those three SNPs, a significant level of the B allele is measured for all three; in this case, the chance that the alleged father is the correct father is low, because the two father haplotypes are A, A, B and B, B, A, while the measured fetal genotype is positive for B at all three SNPs, and the mother could only have contributed an A. If the father genotype was not phased, it would not be possible to rule out this alleged father given this set of measurements. In one embodiment, that determination of the father haplotypes may be determined given the diploid genomic DNA measurements along with haploid genetic measurements made on one or more sperm. The use of more than one sperm can allow the determination of the haplotypes with more accuracy, as well as how many cross overs may have occurred, for each of the chromosomes, along with their locations, during the meiosis that formed the sperm. A method to accomplish this paternal phasing may be found in greater detail in the four Rabinowitz patent applications referenced elsewhere in this document.

In an embodiment, the paternity determination is done exclusively using SNP measurements, and no data from single tandem repeats is used. In an embodiment, the paternity determination is done exclusively using both SNP and STR measurements. The SNP data may be measured using SNP microarrays, or it may be measured by sequencing. The sequencing may be untargeted, or it may be targeted, for example by using circularizing probes that are targeted to a set of polymorphic loci, or it may be use targeted by using capture by hybridization techniques. In some embodiments, the genetic data may be measured by a combination of methods; for example, the parental genetic data may be measured on a SNP microarray while the DNA isolated from maternal serum may be measured using targeted sequencing where capture hybridization probes are used to target the same SNPs as are found on the SNP microarray. In one embodiment a combination of the following types of data may be used to determine whether or not the alleged father is the biological father of the fetus: SNP data, STR data, crossover data, microdeletion data, insertion data, translocation data, or other genetic data.

In an embodiment, the method may comprise the generation of a report disclosing the established paternity of the fetus, or other target individual. In an embodiment, the report may be generated for the purpose of communicating the paternity determination. In an embodiment, the report may comprise a probability that the alleged father is the biological father of the fetus. Some examples of such a report are shown within; FIG. 7 is an example of a report disclosing a paternity exclusion, FIG. 8 is an example of a report disclosing a paternity inclusion and FIG. 9 is an example of a report indicating an indeterminate result. In one embodiment the report may comprise a graph containing a distribution of a paternity related metric for a plurality of unrelated individuals with respect to a given fetus and mother (shown as a grey curve), and an indication of the metric for the alleged father (shown as a triangle). The distribution of unrelated males is different for each case; in these three reports, an actual distribution of the test statistic for the fetus and unrelated males is used here. In an embodiment, the report may also contain an indication that the alleged father is more likely to be part of the distribution of unrelated individuals (e.g. FIG. 7), and therefore the alleged father is established to not be the biological father of the fetus; the fact that the triangle is in the paternity exclusion region of the graph indicates that this is a paternity exclusion. In an embodiment, the report may also contain an indication that the alleged father is more likely to not be part of the distribution of the paternity metric for unrelated individuals (e.g. FIG. 8), and the alleged father is established to be the biological father of the fetus; the fact that the triangle is in the paternity inclusion region of the graph indicates that this is a paternity inclusion. In an embodiment, the report may also contain an indication that the measurements are indeterminate (e.g. FIG. 9); the fact that the triangle is in the “indeterminate result” region of the graph indicates that no conclusion was made with respect to establishing the paternity of the fetus.

In one embodiment of the invention, the determination of whether or not the alleged father is the biological father of the fetus is done without using single tandem repeats (STRs). In one embodiment of the invention, the accuracy of the paternity determination is increased by phasing the parental genotypes. In one embodiment of the invention, the genotypes of one or more of the parents are phased with the use of genetic material from one or more individual related that parent. In one embodiment, the individual related to the parent is the parents father, mother, sibling, son, daughter, brother, sister, aunt, uncle, twin, clone, a gamete from the parent, and combinations thereof.

Another Method for Paternity Determination

In an embodiment, the maternal plasma and optionally the other genetic material may be measured by sequencing, for example using high throughput sequencers such as the HISEQ or MISEQ by ILLUMINA, or the ION TORRENT by LIFE TECHNOLOGIES.

Non-invasive paternity testing can be performed on a maternal blood sample if there is a sufficient concentration of free-floating fetal DNA. In general, the fraction of fetal DNA in most cases, the maternal plasma will be about between 2 percent to 20 percent, though it may be as low as 0.01%, or as high as 40%, partly depending on the gestational age. It has already been demonstrated that this range of fetal fraction is sufficient for paternity testing by a single-hypothesis rejection method using SNP microarrays. High throughput sequencing is a far more precise platform which allows mathematical modeling of the expected measurement response at each SNP, for combinations of mother and child genotypes. In an embodiment, the maternal plasma and optionally the other genetic material may be measured by sequencing, for example using high throughput sequencers such as the HISEQ or MISEQ by ILLUMINA, or the ION TORRENT by LIFE TECHNOLOGIES. Confidences on paternity inclusions or exclusions may then be calculated by using probability and/or estimation theories.

In an embodiment, the method for paternity testing may include the following. For an alleged father, one may calculate the probability of the sequencing data, derived from the plasma, with respect to the two different hypotheses: (1) the alleged father is the correct (biological) father (Hc) and (2) the alleged father is not the correct (biological) father (Hw). The hypothesis that has the higher likelihood or a posteriori is then chosen. In an embodiment, this approach may be combined with a platform model which relates the allele ratio in the plasma to the observed number of sequenced A and B alleles. With the platform model available, it is possible to derive probabilistic likelihoods of the sequenced A and B alleles for each SNP location for each hypothesis.

One complication is that the amount of fetal fraction in the maternal plasma may vary between individuals and over time. In an embodiment, the method may account for this variability. There are several ways to address this type of variability. In an embodiment, the method may explicitly estimate the fetal fraction; in another embodiment, the method may put a prior on the unknown quantity and integrates over all possible values. An embodiment uses a prior that is it as a uniform distribution from 0 to some threshold, e.g. 40%. Any prior may work in theory. An embodiment, calculates likelihoods of various child fractions, either in continuous space or on a finite partition and integrates or sums over the range, respectively.

Consider maternal plasma with fetal fraction, f, and a single SNP where the expected allele ratio present in the plasma is r (based on the maternal and fetal genotypes). In an embodiment, the expected allele ratio is defined as the expected fraction of A alleles in the combined maternal and fetal DNA. For maternal genotype g_(m) and child genotype g_(c), the expected allele ratio is given by equation 1, assuming that the genotypes are represented as allele ratios as well.

r=fg _(c)+(1−f)g _(m)  (1)

The observation at the SNP comprises the number of mapped reads with each allele present, n_(a) and n_(b), which sum to the depth of read d. Assume that quality control measures have been applied to the mapping probabilities such that the mappings and allele observations can be considered correct. A simple model for the observation likelihood is a binomial distribution which assumes that each of the d reads is drawn independently from a large pool that has allele ratio r. Equation 2 describes this model.

$\begin{matrix} {{P\left( {n_{a},{n_{b}❘r}} \right)} = {{p_{bino}\left( {{n_{a};{n_{a} + n_{b}}},r} \right)} = {\begin{pmatrix} {n_{a} + n_{b}} \\ n_{a} \end{pmatrix}{r^{n_{a}}\left( {1 - r} \right)}^{n_{b}}}}} & (2) \end{matrix}$

When the maternal and fetal genotypes are either all A or all B, the expected allele ratio in plasma will be 0 or 1, and p_(bino) will not be well-defined. Additionally, this is not desirable because unexpected alleles are sometimes observed in practice. The binomial model can be extended in a number of ways. In an embodiment, it is possible to use a corrected allele ratio {circumflex over (r)}=1/(n_(a)+n_(b)) to allow a small amount of the unexpected allele to be accounted for. In an embodiment, it is possible to use training data to model the rate of the unexpected allele appearing on each SNP, and use this model to correct the expected allele ratio. When the expected allele ratio is not 0 or 1, the observed allele ratio may not converge to the expected allele ratio due to amplification bias or other phenomena. The allele ratio can then be modeled as a beta distribution centered at the expected allele ratio, leading to a beta-binomial distribution for P(n₈,n_(b)|r) which has higher variance than the binomial.

A general platform model for the response at a single SNP may be defined as F(a, b, g_(c), g_(m), f) (3), or the probability of observing n_(a)=a and n_(b)=b given the maternal and fetal genotypes, which also depends on the fetal fraction through equation 1.

F(a,b,g _(c) ,g _(m) ,f)=P(n _(a) =a,n _(b) =b|g _(c) ,g _(m) ,f)  (3)

Note that it may be feasible to simplify formula (3) by conditioning on a function of g_(c), g_(m) and f e.g. by using r as defined in (1) and the binomial example in (2). The equation for F could then be written

F(a,b,g _(c) ,g _(m) ,f)=P(n _(a) =a,n _(b) =b|g _(c) ,g _(m) ,f)=P(n _(a) =a,n _(b) =b|r(g _(c) ,g _(m) ,f))  (4)

In general the functional form of F may be a binomial distribution, beta-binomial distribution, multivariate Pólya distribution, an empirical distribution estimated from training data, or similar functions as discussed above. In an embodiment, the functional form of F takes different forms depending on the hypothesis for copy number on the chromosome in question.

A Method for the Calculation of the Fetal Fraction

Determining the fraction of fetal DNA that is present in the mixed fraction of DNA may be an integral part of any method for non-invasive prenatal paternity determination, ploidy calling, or allele calling. In some embodiments, the fetal fraction of the mixed sample may be determined using the genotypic data of the mother, the genotypic data of the father, and the measured genotypic data from the mixed sample that contains both maternal and fetal DNA. In the context of paternity testing, and also to a lesser extent in the case of ploidy calling, the identity of the father is not known, and therefore genotypic data of the biological father of a fetus may not be available. In these cases, it is important to have a method for fetal fraction determination that does not require the genotype of the biological father of the fetus. Described herein are several method by which to accomplish the fetal fraction estimate. These methods are described in a general way such that they are appropriate when the genotype of the biological father is available, and when it is not.

For a particular chromosome, suppose we are looking at N SNPs, for which we have the following data:

-   -   A set of NR plasma sequence measurements S=(S₁, . . . ,S_(NR)).         In an embodiment, where we have (A,B) counts for alleles A and B         for each SNP, s can be written as s=((a₁,b₁), . . . ,(a_(N),         b_(N))), where a₁ is the a count on SNP i, b_(i) is the b count         on SNP i, and Σ_(i=1:N)(a_(i)+b_(i))=NR     -   Parent data consisting of:         -   Genotype information: mother G_(m)=(G_(m1), . . . ,G_(mN)),             father G_(f)=(G_(f1), . . . ,G_(fN)), where G_(mi),             G_(fi)∈(AA, AB, BB); and/or         -   Sequence data measurements: NRM mother measurements             s_(m)=(s_(m1), . . . , s_(mnr)), NRF father measurements             S_(f)=(S_(f1), . . . ,S_(fnr)). Similar to above             simplification, if we have (A,B) counts on each SNP             S_(m)=((a_(m1),b_(m1)), . . . ,(a_(mN), b_(mN))),             S_(f)=((a_(f1),b_(f1)), . . . ,(a_(fN), b_(fN)))             Collectively, mother, father child data may be denoted as             D=(G_(m),G_(f),S_(m),S_(f),S). In an embodiment, genotypic             data from both parents are available; in an embodiment,             genotypic data from only the mother is available; in an             embodiment, genotypic data from only the father is             available; in an embodiment, genotypic data from neither             parent is available. In some embodiment, the maternal             genotypic data may be inferred from the genotypic data             measured in the mixed sample. Note that in general, parent             data is desirable and increases the accuracy of the             algorithm, but is not required.

Child fraction estimate {circumflex over (f)} is the expected child fraction given the data:

{circumflex over (f)}=E(cfr|D)=∫f*P(f|D)df

In an embodiment, one may partition the interval of possible child fractions to a set C of finely spaced points and perform the calculations at each point which reduce the above equation to:

$\overset{\hat{}}{f} = {{E\left( f \middle| D \right)} = {\sum\limits_{f \in C}{f*{P\left( f \middle| D \right)}}}}$

P(f|D) is the likelihood of particular child fraction f given the data D. One may further derive using Bayes rule:

P(f|D)−P(D|f)*P(f)

where P(f) is the prior weight of particular child fraction. In an embodiment, this may be derived from uninformed prior(uniform) and may be proportional to the spacing between candidate child fractions in set C.

P(D|cfr) is the likelihood of given data given the particular child fraction, derived under the particular copy number assumptions on the relevant chromosomes. In an embodiment, one may we assume disomy on the chromosomes used. Likelihood of the data on all SNPs is the product of likelihood of data on individual SNPs.

${P\left( D \middle| f \right)} = {{P\left( {\left. D \middle| f \right.,H} \right)} = {\prod\limits_{i}{P\left( {\left. D \middle| f \right.,H,i} \right)}}}$

Where i denotes a particular SNP, for SNP i we have:

${P\left( {\left. D \middle| f \right.,H,i} \right)} = {\sum\limits_{g_{m},g_{f},g_{c}}{{P\left( {\left. D \middle| g_{m} \right.,g_{f},g_{c},f,H,i} \right)}*{P\left( {\left. g_{c} \middle| g_{m} \right.,g_{f},H} \right)}*{P\left( g_{m} \middle| i \right)}*{P\left( g_{f} \middle| i \right)}}}$

where g_(m) are possible true mother genotypes, g_(f) are possible true father genotypes, g_(c) are possible child genotypes, and g_(m), g_(f), g_(c)∈{AA, AB, BB}.

P(g_(m)|i) is the general prior probability of mother genotype g_(m) on SNP i, based on the known population frequency at SNP i, denoted pA_(i). In particular:

p(AA|pA _(i))=(pA _(i))² , p(AB|pA _(i))=2(pA _(i))*(1−pA _(i)), p(BB|pA _(i))=(1−pA _(i))²

Same for p(f|i), father genotype probability.

Let P(g_(c)|g_(m), g_(f), H) denote is the probability of getting true child genotype=c, given parents m, f, and assuming hypothesis H, which we can easily calculate. For example, for a disomy:

parents P(c|m, f, disomy) m f AA AB BB AA AA 1 0 0 AB AA 0.5 0.5 0 BB AA 0 1 0 AA AB 0.5 0.5 0 AB AB 0.25 0.5 0.25 BB AB 0 0.5 0.5 AA BB 0 1 0 AB BB 0 0.5 0.5 BB BB 0 0 1

Let P(D|g_(m), g_(f), g_(c), H, i, f) be the probability of given data D on SNP i, given true mother genotype m, true father genotype f, true child genotype c, hypothesis H for the copy number and child fraction f. It can be broken down into probability of mother, father and child data as follows:

P(D|g _(m) ,g _(f) ,g _(c) ,H,f,i)=P(s _(m) |g _(m) ,i)P(G _(m) |g _(m) ,i)P(s _(f) |g _(f) ,i)P(G _(f) |g _(f) ,i)P(s|g _(m) ,g _(c) ,H,f,i)

The probability of mother illumina genotype data g_(m), at SNP i compared to true genotype g_(m), assuming illumina genotypes are correct, is simply:

${P\left( {\left. G_{m} \middle| g_{m} \right.,i} \right)} = \left\{ \begin{matrix} 1 & {g_{mi} = g_{m}} \\ 0 & {g_{mi} \neq g_{m}} \end{matrix} \right.$

In an embodiment, the probability of mother sequence data at SNP i, in case of counts S_(mi)=(am_(i),bm_(i)), with no extra noise or bias involved, is the binomial probability defined as: P(Sm|,i)=P_(X|m)(am_(i)) where X|m˜Binom(p_(m)(A), am_(i)+bm_(i)) with p_(m)(A) defined as:

m AA AB BB A B nocall p(A) 1 0.5 0 1 0 0.5 A similar equation applies for father probabilities.

Note that it is possible to get an answer without the parent data, especially without the father data. For example if no father genotype data F is available, one can use P(G_(f)|g_(f), i)=1. If no father sequence data S_(f) is available, one can use P(S_(f)|g_(f),i)=1. In an embodiment, information from different chromosomes is aggregated using averages, weighted average or a similar function.

Another Method for the Calculation of the Fetal Fraction

Another method for determining the fraction of fetal DNA in a mixture of DNA is described here. In one embodiment, a version of a maximum likelihood estimate of the fetal fraction f for a paternity test, ploidy test, or other purpose, may be derived without the use of paternal information. Define S₀ as the set of SNPs with maternal genotype 0 (AA), S_(0.5) as the set of SNPs with maternal genotype 0.5 (AB) and S₁ as the set of SNPs with maternal genotype 1 (BB). The possible fetal genotypes on S₀ are 0 and 0.5, resulting in a set of possible allele ratios R₀(f)={0,f/2}. Similarly, R_(0.5)={0.5−f, 0.5, 0.5+f} and R₁(f)={1−f/2, 1}. All or any subset of the sets S₀, S_(0.5) and S₁ can be used to derive a child fraction estimate.

Define N_(a0) and N_(b0) as the vectors formed by sequence counts for SNPs in S₀, N_(0.5) and Nb_(0.5) similarly for S_(0.5), and N_(a1) and N_(b1) similarly for S₁. The maximum likelihood estimate {circumflex over (f)} of f, using all maternal genotype sets, is defined by equation 4.

{circumflex over (f)}=arg max_(f) P(N _(a0) ,N _(b0) |f)P(Na _(0.5) ,Nb _(0.5) |f)P(N _(a1) ,N _(b1) |f)  (4)

Assuming that the allele counts at each SNP are independent conditioned on the SNP's plasma allele ratio, the probabilities can be expressed as products over the SNPs in each set:

P(N _(a0) ,N _(b0) |f)=ΠsϵS ₀ P(n _(as) ,n _(bs) |f)

P(N _(a1) ,N _(b1) |f)=ΠsϵS ₁ P(n _(as) ,n _(bs) |f)

where n_(as), n_(bs) are the counts on SNPs s.

The dependence on f is through the sets of possible allele ratios R₀(f), R_(0.5)(f) and R₁(f). The SNP probability P(n_(as), n_(bs)|f) can be approximated by assuming the maximum likelihood genotype conditioned on f. At reasonably high fetal fraction and depth of read, the selection of the maximum likelihood genotype will be high confidence. For example, at fetal fraction of 10 percent and depth of read of 1,000, consider a SNP where the mother has genotype 0. The expected allele ratios are 0 and 5 percent, which will be easily distinguishable at sufficiently high depth of read. Substitution of the estimated child genotype into equation 5 results in the complete equation (6) for the fetal fraction estimate.

$\overset{\hat{}}{f} = {\arg\mspace{14mu}{\max_{f}\begin{bmatrix} {{\Pi_{{s\epsilon S}_{0}}\left( {\max\limits_{r_{s}{{\epsilon R}_{0}{(f)}}}{P\left( {n_{as},{n_{bs}❘r_{s}}} \right)}} \right)}{\Pi_{{s\epsilon S}_{0.5}}\left( {\max\limits_{r_{s}{{\epsilon R}_{0.5}{(f)}}}{P\left( {n_{as},{n_{bs}❘r_{s}}} \right)}} \right)}} \\ {\Pi_{{s\epsilon S}_{1}}\left( {\max\limits_{r_{s}{{\epsilon R}_{1}{(f)}}}{P\left( {n_{as},{n_{bs}❘r_{s}}} \right)}} \right.} \end{bmatrix}}}$

The fetal fraction must be in the range [0, 1] and so the optimization can be easily implemented by a constrained one-dimensional search.

Another method would be to sum over the possible genotypes at each SNP, resulting in the following expression (7) for P(n_(a), n_(b)|f) for a SNP in S₀. The prior probability P(r) could be assumed uniform over R₀(f), or could be based on population frequencies. The extension to groups S_(0.5) and S₁ is trivial.

P(n _(a) ,n _(b) |f)=Σ_(r∈R) ₀ _((f)) P(n _(a) ,n _(b) |r)P(r)  (7)

Derivation of Probabilities

A confidence can be calculated from the data likelihoods of the two hypotheses H_(tf) i.e. the alleged father is the biological father and H_(wf) i.e. the alleged father is not the biological father. The objective is to calculate P(H|D) i.e. probability of hypothesis given data, for each hypothesis and infer which hypothesis is more likely. In one embodiment, this may be done using Bayes rule: P(H|D)˜P(D|H)*P(H) where P(H) is the prior weight of the hypothesis, and where P(D|H) is the likelihood of data given the hypothesis.

Consider P(D|H,f) i.e. the likelihood of data given hypothesis for a particular child fraction. If a distribution on child fraction is available, it is possible to derive

P(D|H)=∫P(D,f|H)df

and further,

P(D|H)=∫P(D|H,f)P(f|H)df

Note that P(f|H) is independent of the hypothesis i.e. P(f|H)=P(f) since the child fraction is the same regardless of whether the alleged father is the biological father or not, and any reasonable prior P(f) could be chosen e.g. a uniform prior from 0 to 50% child fraction. In an embodiment, it is possible to use only one child fraction, {circumflex over (f)}. In this case,

P(D|H)=P(D|H,{circumflex over (f)})

Consider the likelihood P(D|H,f). The likelihood of each hypothesis is derived based on the response model, the estimated fetal fraction, the mother genotypes, the alleged father genotypes, allele population frequencies, the plasma allele counts and SNPs. Let D represent the data as defined before.

In an embodiment, it is assumed that the observation at each SNP is independent conditioned on the plasma allele ratio, thus the likelihood of a paternity hypothesis is the product of the likelihoods on the SNPs:

${P\left( {\left. D \middle| H \right.,f} \right)} = {\prod\limits_{SNPsi}{P\left( {\left. D \middle| H \right.,f,i} \right)}}$

The following equations describe how one may derive the likelihood for a single SNP i and a single child fraction f. Equation 8 is a general expression for the likelihood of any hypothesis H, which will then be broken down into the specific cases of H_(tf) and H_(wf). Note that genotypes, g_(m), g_(tf), g_(df), and g_(c), take values in {AA,AB,BB} which translates to {0,0.5,1} where AA=0, AB=0.5, BB=1. Also, g_(tf) denote the genotypes of the true father and g_(df) denote the genotypes represented by the data provided for the fater. In case of H_(tf), g_(tf) and g_(df) are equivalent.

P(D|f,H,i)=Σ_(g) _(m) _(,g) _(tf) _(,g) _(df) _(,g) _(c) _(∈{0,0.5,1}) P(D|g _(m) ,g _(df) ,g _(c) ,f,H,i)*P(g _(c) |g _(m) ,g _(tf) ,H)*P(g _(ml) |i)*P(g _(tf) |i)*P(g _(df) |i)  (8)

In the case of the hypothesis H_(tf), the alleged father is the biological father and the fetal genotypes are inherited from the maternal genotypes and alleged father genotypes. The equation above simplifies to:

${P\left( {{D❘f},{H = H_{tf}},i} \right)} = {\underset{g_{m},g_{tf},{g_{c}\epsilon{\{{0,0.5,1}\}}}}{\Sigma}{P\left( {{D❘g_{m}},g_{tf},g_{c},f,{H = H_{tf}},i} \right)}*{P\left( {{g_{c}❘g_{m}},g_{tf},{H = H_{tf}}} \right)}*{P\left( {g_{m}❘i} \right)}*{P\left( {g_{tf}❘i} \right)}}$

Further,

P(g _(c) |g _(m) ,g _(tf) ,H=H _(tf))=P(g _(c) |g _(m) ,g _(tf))

and

P(D|g _(m) ,g _(tf) ,g _(c) ,f,H=H _(tf) ,i)=P(s _(m) |g _(m) ,i)P(G _(m) |g _(m) ,i)P(s _(f) |g _(tf) ,i)P(G _(f) |g _(tf) ,i)P(s|g _(m) ,g _(c) ,f,i)  (9)

In the case of H_(wf), the alleged father is not the biological father. One estimate of the true father genotypes may be generated using the population frequencies at each SNP. Thus, the probabilities of child genotypes may be determined by the known mother genotypes and the population frequencies, i.e. the data do not provide additional information on the genotypes of the biological father. In this case, the equation above does not further simplify and stays as:

${P\left( {\left. D \middle| f \right.,{H = H_{wf}},i} \right)} = {\sum\limits_{g_{m},g_{tf},g_{df},{g_{c}\epsilon{\{{0,{0.5},1}\}}}}{{P\left( {\left. D \middle| g_{m} \right.,g_{df},g_{c},f,{H = H_{wf}},i} \right)}*{P\left( {\left. g_{c} \middle| g_{m} \right.,g_{tf},{H = H_{wf}}} \right)}*{P\left( g_{m} \middle| i \right)}*{P\left( g_{tf} \middle| i \right)}*{P\left( g_{df} \middle| i \right)}}}$

Further,

P(g _(c) |g _(m) ,g _(tf) ,H=H _(wf))=P(g _(c) |g _(m) ,g _(tf))

where the only information on g_(tf) are the population priors and:

P(D|g _(m) ,g _(tf) ,g _(c) ,f,H=H _(wf) ,i)=P(s _(m) |g _(m) ,i)P(G _(m) |g _(m) ,i)P(s _(f) |g _(df) ,i)P(G _(f) |g _(df) ,i)P(s|g _(m) ,g _(c) ,f,i)  (10)

In both expressions of the likelihoods, P(D|f,H,i), i.e. for both hypotheses, the response model, P(s|g_(m),g_(df),g_(c),f,H) is generalized. Specific examples are mentioned elsewhere in the document in discussions on general platform models. Some examples for the response model include the binomial distribution, beta-binomial distribution, multivariate Pólya distribution, an empirical distribution estimated from training data, or similar functions as discussed above.

In some embodiments, the confidence C, on correct paternity can be calculated from the likelihoods P(D|H_(tf)) and P(D|H_(wf)). In an embodiment this calculation may be calculated using Bayes rule as follows:

$C_{p} = \frac{{P\left( D \middle| H_{tf} \right)}{P\left( H_{tf} \right)}}{{{P\left( D \middle| H_{tf} \right)}{P\left( H_{tf} \right)}} + {{P\left( D \middle| H_{wf} \right)}{P\left( H_{wf} \right)}}}$

or written for a more specific case as a product over SNPs of the two likelihoods:

$C_{p} = \frac{\Pi_{i}{P\left( {\left. D \middle| H_{tf} \right.,G_{ms},G_{tf},f} \right)}}{{\Pi_{s}{P\left( {n_{as},{n_{bs}❘H_{t}},G_{ms},G_{tf},f} \right)}} + {\Pi_{s}{P\left( {n_{as},{n_{bs}❘H_{f}},G_{ms},G_{tf},f} \right)}}}$

In another embodiment the confidence may be calculated as follows:

$C_{p} = \frac{P\left( D \middle| H_{tf} \right)}{{P\left( D \middle| H_{tf} \right)} + {P\left( D \middle| H_{wf} \right)}}$

Other reasonable functions of the likelihoods are also possible.

EXPERIMENTAL SECTION Experiment 1

Twenty one pregnant women with confirmed paternity and gestational ages between 6 and 21 weeks were enrolled. Participants voluntarily donated blood as part of our IRB approved research program, and were drawn from IVF centers, OB offices, and the general population in different locations in the U.S. Cell free DNA (ffDNA) isolated from maternal plasma, along with DNA from the mother and alleged father, were amplified and measured using a SNP array. An informatics method disclosed herein was used to exclude or include paternity for 21 correct fathers and 36,400 incorrect fathers by comparing each alleged father against a reference distribution generated from a set of over 5,000 unrelated individuals. 20 out of 21 samples had sufficient fetal DNA to return results. Twenty of twenty (100%) of paternity inclusions were correct. 36,382 of 36,382 paternity exclusions were correct (100%), with 18 “no calls” due to intermediate genetic similarity. There were no miscalls.

The population was made up of couples who donated their blood for prenatal research. The women had to have singleton pregnancies, be in the first or second trimester, and have confirmed paternity. Blood samples were collected from women using CELL-FREE blood tubes (STRECK) containing white blood cell preservative, and genetic samples were collected from the father, either as a blood (EDTA) or buccal sample. Written informed consent was obtained from all participants, and the genetic samples were collected from patients enrolled in an IRB approved study.

Mother blood was centrifuged to isolate the buffy coat and the serum. The genomic DNA in the maternal and alleged paternal buffy coat and the DNA in the maternal serum were prepared for analysis and run on ILLUMINA INFINIUM CYTO12 SNP arrays using standard protocols. Briefly, serum DNA was isolated using QIAGEN CIRCULATING NUCLEIC ACID kit and eluted in 45 ul buffer according to manufacturer's instructions. Twenty microliter eluate was used in a blunt ending reaction in 1×NEB 4 buffer, 0.42 mM dNTP and 2.5 U T4 DNA Polymerase (NEW ENGLAND BIOLABS), incubated at 20 C for 30 min, then 75 C for 15 min. Three microliter ligation mixture (0.5 ul 10×NEB 4, 1 ul 10 mM ATP, 1 ul T4 PNK (NEW ENGLAND BIOLABS), 0.5 ul T4 DNA Ligase (NEW ENGLAND BIOLABS)) was added and samples incubated at 16 C for 24 hours, then 75 C for 15 min. The sample was transferred to the standard ILLUMINA INFINIUM assay along with the maternal and alleged paternal genomic DNA. In short, 24 ul DNA was whole genome amplified at 37 C for 20-24 hours followed by fragmentation and precipitation. The precipitate was then resuspended in hybridization buffer, heat denatured and transferred to Cyto12 SNP arrays using a TECAN EVO. The arrays were incubated at 48 C for at least 16 hours, X-Stained (INFINIUM II CHEMISTRY) and washed in the TECAN EVO, and finally scanned. Array intensities were extracted using BEADSTIDUO (ILLUMINA).

The disclosed informatics method generated a test statistic that measured the degree of genetic similarity between the fetus and another individual. This test statistic was calculated for both the alleged father and a set of over 5,000 unrelated individuals. A single hypothesis rejection test then determined whether the statistic calculated for the alleged father could be excluded from the distribution formed by the unrelated reference individuals. If the alleged father could be rejected from the unrelated set, then a paternity inclusion resulted; otherwise, paternity was excluded. For the 20 samples with sufficient DNA, the paternity test was run against 20 correct fathers, and for 1,820 randomly selected incorrect fathers.

A paternity inclusion was called when the p-value of the alleged father's test statistic on the distribution of unrelated individuals was less than 104. This means that, in theory, no more than one out of 10,000 unrelated individuals are expected to show as much genetic similarity to the fetus. A “no call” was called where the p-value is between 10 and 0.02. An “insufficient fetal DNA” call was made when the fetal DNA made up less than 2% of the plasma DNA. The set of unrelated individuals used to generate the expected distribution was composed of individuals from a wide variety of racial backgrounds, and the paternity inclusion or exclusion determination was recalculated for sets of unrelated individuals of different races, including the race indicated for the alleged father. The inclusion and exclusion results were automatically generated by the algorithm, and no human intervention was necessary.

In conclusion, twenty one samples of maternal blood with known paternity were tested. Twenty out of 21 samples returned results, while one had insufficient fetal DNA for analysis; this sample was drawn from a woman at 8 weeks gestational age. Twenty of twenty (100%) results had the correct paternity confirmed, each with a p-value of <104. Each of the 20 samples with sufficient fetal fraction was tested against a random set of 1,820 incorrect fathers, for a total of 36,400 individual paternity tests. 36,382 of these analyses returned a result; 36,382 of 36,382 (100%) correctly had the paternity excluded with a p-value of greater than 104, and 18 of 36,400 (0.05%) were called “no call”, with a p-value between 104 and 0.02. There were no incorrect paternity exclusions or inclusions.

Nine of 21 samples had confirmed paternity due to control of fertilization during IVF with correct paternity confirmed after fertilization through pre-implantation genetic diagnosis. Twelve samples had paternity confirmed by independent paternity testing of fetal/child genomic DNA, conducted by DNA Diagnostic Center, Fairfield, Ohio.

Experiment 2

In one experiment, four maternal plasma samples were prepared and amplified using a hemi-nested 9,600-plex protocol. The samples were prepared in the following way: between 15 and 40 mL of maternal blood were centrifuged to isolate the buffy coat and the plasma. The genomic DNA in the maternal and was prepared from the buffy coat and paternal DNA was prepared from a blood sample or saliva sample. Cell-free DNA in the maternal plasma was isolated using the QIAGEN CIRCULATING NUCLEIC ACID kit and eluted in 45 uL TE buffer according to manufacturer's instructions. Universal ligation adapters were appended to the end of each molecule of 35 uL of purified plasma DNA and libraries were amplified for 7 cycles using adaptor specific primers. Libraries were purified with AGENCOURT AMPURE beads and eluted in 50 ul water.

3 ul of the DNA was amplified with 15 cycles of STA (95° C. for 10 min for initial polymerase activation, then 15 cycles of 95° C. for 30 s; 72° C. for 10 s; 65° C. for 1 min; 60° C. for 8 min; 65° C. for 3 min and 72° C. for 30 s; and a final extension at 72° C. for 2 min) using 14.5 nM primer concentration of 9600 target-specific tagged reverse primers and one library adaptor specific forward primer at 500 nM.

The hemi-nested PCR protocol involved a second amplification of a dilution of the first STAs product for 15 cycles of STA (95° C. for 10 min for initial polymerase activation, then 15 cycles of 95° C. for 30 s; 65° C. for 1 min; 60° C. for 5 min; 65° C. for 5 min and 72° C. for 30 s; and a final extension at 72° C. for 2 min. using reverse tag concentration of 1000 nM, and a concentration of 16.6 u nM for each of 9600 target-specific forward primers.

An aliquot of the STA products was then amplified by standard PCR for 10 cycles with 1 uM of tag-specific forward and barcoded reverse primers to generate barcoded sequencing libraries. An aliquot of each library was mixed with libraries of different barcodes and purified using a spin column.

In this way, 9,600 primers were used in the single-well reactions; the primers were designed to target SNPs found on chromosomes 1, 2, 13, 18, 21, X and Y. The amplicons were then sequenced using an ILLUMINA GAIIX sequencer. Per sample, approximately 3.9 million reads were generated by the sequencer, with 3.7 million reads mapping to the genome (94%), and of those, 2.9 million reads (74%) mapped to targeted SNPs with an average depth of read of 344 and a median depth of read of 255. The fetal fraction for the four samples was found to be 9.9%, 18.9%, 16.3%, and 21.2%

Relevant maternal and paternal genomic DNA samples amplified using a semi-nested 9600-plex protocol and sequenced. The semi-nested protocol is different in that it applies 9,600 outer forward primers and tagged reverse primers at 7.3 nM in the first STA. Thermocycling conditions and composition of the second STA, and the barcoding PCR were the same as for the hemi-nested protocol.

The sequencing data was analyzed using informatics methods disclosed herein and each of a set of ten unrelated males from a reference set were determined to not be the biological father of each of the gestating fetuses.

Experiment 3

In one experiment 45 sets of cells were amplified using a 1,200-plex semi-nested protocol, sequenced, and ploidy determinations were made at three chromosomes. Note that this experiment is meant to simulate the conditions of performing paternity testing on single fetal cells obtained from maternal blood, or on forensic samples where a small amount of DNA from the child is present. 15 individual single cells and 30 sets of three cells were placed in 45 individual reaction tubes for a total of 45 reactions where each reaction contained cells from only one cell line, but the different reactions contained cells from different cell lines. The cells were prepared into 5 ul washing buffer and lysed the by adding 5 ul ARCTURUS PICOPURE lysis buffer (APPLIED BIOSYSTEMS) and incubating at 56° C. for 20 min, 95° C. for 10 min.

The DNA of the single/three cells was amplified with 25 cycles of STA (95° C. for 10 min for initial polymerase activation, then 25 cycles of 95° C. for 30 s; 72° C. for 10 s; 65° C. for 1 min; 60° C. for 8 min; 65° C. for 3 min and 72° C. for 30 s; and a final extension at 72° C. for 2 min) using 50 nM primer concentration of 1200 target-specific forward and tagged reverse primers.

The semi-nested PCR protocol involved three parallel second amplification of a dilution of the first STAs product for 20 cycles of STA (95° C. for 10 min for initial polymerase activation, then 15 cycles of 95° C. for 30 s; 65° C. for 1 min; 60° C. for 5 min; 65° C. for 5 min and 72° C. for 30 s; and a final extension at 72° C. for 2 min) using reverse tag specific primer concentration of 1000 nM, and a concentration of 60 nM for each of 400 target-specific nested forward primers. In the three parallel 400-plex reactions the total of 1200 targets amplified in the first STA were thus amplified.

An aliquot of the STA products was then amplified by standard PCR for 15 cycles with 1 uM of tag-specific forward and barcoded reverse primers to generate barcoded sequencing libraries. An aliquot of each library was mixed with libraries of different barcodes and purified using a spin column.

In this way, 1,200 primers were used in the single cell reactions; the primers were designed to target SNPs found on chromosomes 1, 21 and X. The amplicons were then sequenced using an ILLUMINA GAIIX sequencer. Per sample, approximately 3.9 million reads were generated by the sequencer, with 500,000 to 800,000 million reads mapping to the genome (74% to 94% of all reads per sample).

Relevant maternal and paternal genomic DNA samples from cell lines were analyzed using the same semi-nested 1200-plex assay pool with a similar protocol with fewer cycles and 1200-plex second STA, and sequenced.

The sequencing data was analyzed using informatics methods disclosed herein and each of a set of ten unrelated males from a reference set were determined to not be the biological father of the target individual for each of the 45 cells.

DNA from Children from Previous Pregnancies in Maternal Blood

One difficulty to non-invasive prenatal paternity testing is differentiating fetal cells from the current pregnancy from fetal cells from previous pregnancies. Some believe that genetic matter from prior pregnancies will go away after some time, but conclusive evidence has not been shown. In an embodiment of the present disclosure, it is possible to determine fetal DNA present in the maternal blood of paternal origin (that is, DNA that the fetus inherited from the father) using the PARENTAL SUPPORT™ (PS) method, and the knowledge of the paternal genome. This method may utilize phased parental genetic information. It is possible to phase the parental genotype from unphased genotypic information using grandparental genetic data (such as measured genetic data from a sperm from the grandfather), or genetic data from other born children, or a sample of a miscarriage. One could also phase unphased genetic information by way of a HapMap-based phasing, or a haplotyping of paternal cells. Successful haplotyping has been demonstrated by arresting cells at phase of mitosis when chromosomes are tight bundles and using microfluidics to put separate chromosomes in separate wells. In another embodiment it is possible to use the phased parental haplotypic data to detect the presence of more than one homolog from the father, implying that the genetic material from more than one child is present in the blood. By focusing on chromosomes that are expected to be euploid in a fetus, one could rule out the possibility that the fetus was afflicted with a trisomy. Also, it is possible to determine if the fetal DNA is not from the current father, in which case one could use other methods such as the triple test to predict genetic abnormalities.

There may be other sources of fetal genetic material available via methods other than α blood draw. In the case of the fetal genetic material available in maternal blood, there are two main categories: (1) whole fetal cells, for example, nucleated fetal red blood cells or erythroblats, and (2) free floating fetal DNA. In the case of whole fetal cells, there is some evidence that fetal cells can persist in maternal blood for an extended period of time such that it is possible to isolate a cell from a pregnant woman that contains the DNA from a child or fetus from a prior pregnancy. There is also evidence that the free floating fetal DNA is cleared from the system in a matter of weeks. One challenge is how to determine the identity of the individual whose genetic material is contained in the cell, namely to ensure that the measured genetic material is not from a fetus from a prior pregnancy. In an embodiment of the present disclosure, the knowledge of the maternal genetic material can be used to ensure that the genetic material in question is not maternal genetic material. There are a number of methods to accomplish this end, including informatics based methods such as PARENTAL SUPPORT™, as described in this document or any of the patents referenced in this document.

In an embodiment of the present disclosure, the blood drawn from the pregnant mother may be separated into a fraction comprising free floating fetal DNA, and a fraction comprising nucleated red blood cells. The free floating DNA may optionally be enriched, and the genotypic information of the DNA may be measured. From the measured genotypic information from the free floating DNA, the knowledge of the maternal genotype may be used to determine aspects of the fetal genotype. These aspects may refer to ploidy state, and/or a set of allele identities. Then, individual nucleated cells that are presumably or possible fetal in origin may be genotyped using methods described elsewhere in this document, and other referent patents, especially those mentioned in this document. The knowledge of the maternal genome would allow one to determine whether or not any given single blood cell is genetically maternal. And the aspects of the fetal genotype that were determined as described above would allow one to determine if the single blood cell is genetically derived from the fetus that is currently gestating. In essence, this aspect of the present disclosure allows one to use the genetic knowledge of the mother, and possibly the genetic information from other related individuals, such as the father, along with the measured genetic information from the free floating DNA found in maternal blood to determine whether an isolated nucleated cell found in maternal blood is either (a) genetically maternal, (b) genetically from the fetus currently gestating, or (c) genetically from a fetus from a prior pregnancy.

All patents, patent applications, and published references cited herein are hereby incorporated by reference in their entirety. It will be appreciated that several of the above-disclosed and other features and functions, or alternatives thereof, may be desirably combined into many other different systems or applications. Various presently unforeseen or unanticipated alternatives, modifications, variations, or improvements therein may be subsequently made by those skilled in the art which are also intended to be encompassed by the following claims. 

What is claimed is:
 1. A method for nested PCR amplification, comprising: isolating cell-free DNA from a biological sample and ligating adaptors to the isolated cell-free DNA, wherein the adaptors each comprise a universal priming site; performing a first PCR to simultaneously amplify at least 10 target loci present in the cell-free DNA using a universal primer that binds the universal priming site, and at least 10 target-specific primers; and performing a second, nested PCR to amplify the at least 10 target loci using the universal primer and at least 10 inner target-specific primers to obtain amplified DNA, wherein primer binding sites of the inner target-specific primers of the second PCR are internal to primer binding sites of the target-specific primers of the first PCR.
 2. The method of claim 1, wherein the biological sample is a blood, plasma, serum, or urine sample.
 3. The method of claim 1, wherein the biological sample is a plasma sample.
 4. The method of claim 1, further comprising sequencing the amplified DNA.
 5. The method of claim 1, wherein the adaptors are looped ligation adaptors.
 6. The method of claim 1, wherein the at least 10 inner target-specific primers comprise a target-specific region that binds to the target loci and a priming site 5′ to the target-specific region that does not bind to the target loci.
 7. The method of claim 6, further comprising performing a third PCR using a forward primer that comprises the sequence of the universal primer and a reverse primer that binds to the priming site on the amplified DNA to generate sample barcoded DNA, wherein the forward primer further comprises a sample barcode sequence.
 8. The method of claim 7, further comprising sequencing the sample barcoded DNA.
 9. The method of claim 8, wherein the sequencing is multiplex sequencing of multiple samples.
 10. The method of claim 6, wherein at least one of the primers comprises a partial sequencing tag.
 11. The method of claim 1, wherein the first PCR comprises simultaneously amplifying at least 50 target loci using the universal primer and at least 50 target-specific primers in a first reaction volume.
 12. The method of claim 11, wherein the second, nested PCR comprises amplifying the at least 50 target loci using the universal primer and at least 50 target-specific primers in a plurality of sub-reactions.
 13. The method of claim 1, wherein the first PCR comprises simultaneously amplifying at least 50 target loci using the universal primer and at least 50 target-specific primers in a single reaction volume.
 14. The method of claim 1, wherein the first PCR comprises simultaneously amplifying at least 100 target loci using the universal primer and at least 100 target-specific primers in a single reaction volume.
 15. The method of claim 1, wherein the concentration of each target-specific primer of the first PCR and/or second PCR is less than 20 nM.
 16. The method of claim 1, wherein the concentration of each target-specific primer of the first PCR and/or second PCR is less than 10 nM.
 17. The method of claim 1, wherein the target loci are SNP loci.
 18. The method of claim 1, wherein the cell-free DNA comprises DNA from a tumor.
 19. The method of claim 1, wherein the cell-free DNA comprises DNA from a transplant.
 20. The method of claim 1, wherein the adaptors each comprise a molecular barcode.
 21. The method of claim 1, wherein the isolated cell-free DNA is pre-amplified before ligating the adaptors.
 22. The method of claim 21, wherein the pre-amplification of the isolated cell-free DNA is a linear amplification.
 23. The method of claim 22, wherein the linear amplification is multiplexed linear amplification.
 24. The method of claim 1, wherein the at least 10 target-specific primers are selected from a plurality of primers, wherein the at least 10 target-specific primers are selected by identifying primers in the plurality of primers that are likely to cause unproductive side reactions, and not including the identified primers in the at least 10 target-specific primers.
 25. The method of claim 24, wherein sequencing of DNA amplified using the plurality of primers is performed to generate sequencing data, and wherein the at least 10 target-specific primers are selected by identifying the primers in the plurality of primers with the most non-mapping primer dimer side reactions, and not including the identified primers in the at least 10 target-specific primers.
 26. The method of claim 24, wherein the binding energies of various primer combinations are calculated, and wherein the at least 10 target-specific primers are selected by identifying the primers in the plurality of primers with the highest binding energies, and not including the identified primers in the at least 10 target-specific primers.
 27. The method of claim 1, wherein the at least 10 inner target-specific primers are selected from a plurality of primers, wherein the at least 10 inner target-specific primers are selected by identifying primers in the plurality of primers that are likely to cause unproductive side reactions, and not including the identified primers in the at least 10 inner target-specific primers.
 28. The method of claim 27, wherein sequencing of DNA amplified using the plurality of primers is performed to generate sequencing data, and wherein the at least 10 inner target-specific primers are selected by identifying the primers in the plurality of primers with the most non-mapping primer dimer side reactions, and not including the identified primers in the at least 10 inner target-specific primers.
 29. The method of claim 27, wherein the binding energies of various primer combinations are calculated, and wherein the at least 10 inner target-specific primers are selected by identifying the primers in the plurality of primers with the highest binding energies, and not including the identified primers in the at least 10 inner target-specific primers.
 30. The method of claim 1, wherein at least 80% of the amplified DNA maps to the target loci. 